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Addendum to a Guide to Tables on 
Punched Cards 


A Guide to Tables on Punched Cards was published in October, 1951 
[MTAC, v. 5, p. 185-212]. In July, 1952 [MTAC, v. 6, p. 204-205], an 
announcement was made about a forthcoming addendum to this guide, and 
information was requested concerning any new contributions or necessary 
corrections. This addendum and list of corrections is now given below. In 
; addition, information is supplied in Section 2 about keypunching of im- 
portant tables that is now in progress at various laboratories. 

It might be interesting to observe that the supply of new tables on 
punched cards is rather small, compared with the activity in mathematical 
computing that has been stimulated by the increasing number of high-speed 
calculators. It is possible that this foreshadows the replacement of punched 
cards by more satisfactory tools, which may become available even to labora- 


a tories of relatively small size. For the present, however, punched cards are 
) still very much in use and it is deemed desirable to keep the record of them 
up-to-date. 
Additions: 
Source Available at Description of Tables 
[(18)] (18)C 2.6 vx, Vx/10: 
x = 0(.01)6(.1)14.6; 5D 
[(18)] (18)C 2.6 Vx, ¥x/10: x = 0(.001)2; 5D 
[(18)] (18)C 2.6 Vx:x = 1(1)9999; 9D; 4000 cards 
(x = 0001, 1001, 2001 on same 
card, etc.) 
[(24)] (24)C 2.9 Sums related to powers of complex 
numbers. 


U, = 3"/n!; T, = U,/n 


n=l 


Pi t+ iQe = 

n=l 
Ri, Sz, Pe, Qe, R = 1(1)4; 10D 
2=x+ ty;x,y = 0(.1)3.1; 
Useful for computing the expo- 


nential functions and their inte- 
grals in the complex plane. 


Available at 


(24)C, 
(73)UMT 145[D,F] 


Source 


[(24)] 


[(67)] (67) 
[(67)] (67) 
[(67)] (67) 
[(67)] (67) 
[(67)] (67) 
[(67)] (67) 


[(23)-(18)] (18)C 
[(18)] (18)C 


[(16)-(4)] (4) 


[(25)] (25) 


C(5)] 


[(66)-(4)] (4) 
[(66)-(4)] (4) 


[(66)-(4)]@) { 


(18)C 


7.1 


7.4 


7.4 


7.5 


7.5 


9.10 


9.15 


10.0 
10.0 


14.61 


14.9 


17.1 


18.1 
18.3 


18.3 
18.4 
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Description of Tables 


Cyclotomic cosines 


2 cos (24k/p): 

iw 10( 25+ 20D 
p an odd prime < 100. 
517 values 


(Computed under supervision of 
D. H. LEHMER.) 
sin x, cos x, A, A/60: 
x = 0(1')90°; 8D 
(See Note 1 at end of section. ) 
sin x, cos x, D!, D?: 
x = 0’’(1000’)129600” 8D 
(See Note 1 at end of section.) 
sin x, cos x, A: 
x = 0*(100*)172800"; 5D 
sin x, cos x, A: 
x = 0(0*.01)11*.99; 4D 
= arc sin x, D'!, D?: 
x = 0(.01)0.52; y to 6D of a degree 
= arc tan x: 
x = 0(.01)1; y to 7D of an hour 
e-*: x = 0(.01)30; 14D 
e-*: x = 0(.001)2; 
also x = 0(.01)6(.1)14.6; 5D 

1 
* =U+iV:2=x + ty; 
x = — .5(.01).5; y = 0(.01)1; 6D 
Incomplete Beta function 
B-(p,g):x = 1,p < = 0.05 
(.05)10; 6S 
x = $;p = g = .05(.05)10; 6S 
See Bureau of Mines R I 4917, in 
press, for details and tabular 
entries. 
Jo(x), Ji(x), A, A’, A’, At: 
x = 25.01(.01)99.9; 18D 
Io(x), 6: x = 0(.001)5; 9S 
Ko(x), &: x = .01(.01)5; 8S 
Ip(x), &, e* K(x), 
x = 5(.01)10(.1)20; 8S 
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Source Available at 
[(24) ] (24)C 


[(24)—-(12)] (24), (12)C 


[(4)] (4)C 


[(1S)] (15) 


[(18)] (18)C 


[(4)-(4)] (4) 


[(14)-(4)] (4) 


Description of Tables 


20.294 f(x) = f Ai(—t)dt 


22.56 


22.68 


25.1 


25.1 


25.21 


25.31 


Fe) = 


where Ai(#) is the familiar Airy 
integral: 

x = — 2(.01)5; 8D in f(x); 

7D and 8D in F(x), with the 
eighth place uncertain. Computed 
by E. Ossorne. 


e-* L,(2x), where L,(2x) is the 
Laguerre polynomial : 

n = 0(1)12, 

x = 0(.1).5(.5)10(1)20; 10D 

nm = 13(1)20; x = 0(.5)10(1)20; 
10D 

Computed under supervision of 
C. Lanczos. 


Fermi functions (see definitions 
in note 2 at end of section). 

Table 1. 

do(p, 2): |z| = 2(1)96; 

p = 0(.1)1(.2)3(.5)5(1)9(2)15 

|z| = 2(1)25; p = 20, 25; to 
about 4S 

Table 2. 

F, 


WF: v = 1, 2, 3; same range as 


Table 1, except that p > 0.1. 


Random digits. Random with 
respect to card column, as well as 
with respect to digits in col. Made 
by JACK SHERMAN. See Note 4. 


8-digit random numbers com- 
puted by Lehmer’s rule—9 num- 
bers per card, 3,000 cards 


Random variates from distribu- 
tion exp(—y). [A number r was 
chosen at random from 4-decimal 
numbers. Then — logr = y]. 
140,000 cards. 


Random direction cosines from a 
spherically symmetric distribu- 
tion. Three-dimensional, /, m, n, 
to 3 decimals. 70,000. 


| 
| 
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Source Available at Description of Tables 

[(25)] [25] 26.1 Tables relating to hydrodyna- 
mics 
Tables of velocity of steady lami- 
nar flow in channels of rectangular 
cross-section (see definition in 
note 3 at end of section). 
w(x, y, ¢) to 6D for x, y, o ranging 
between 0 and 1. 

Northrop (73)UMT138[F] 27.1 Products of powers of small primes 

N = 2¢ 3°5¢74:a4 = O(1)11; 

b = 0(1)8; c = 0(1)5; d = 0(1)4; 
exact. 
Computed by R. A. JoHNson, 
Northrop Aircraft, Inc., Haw- 
thorne, Calif. 


(24) ] (24)C, 27.7 Kloosterman sums 
UMT 146[F,L] p-1 
Sp(k) = Lexp{2xi(kn + n)/p} 


nt = 1(mod p) 


k=1(1) (p-—1), p a prime 
< 100, 19D. 

1034 cards. Computed under 
supervision of D. H. LEHMER. 


[(18)] (18)C 28.1 Astronomy. Julian date to calen- 
dar date 1639 (30 Julian days) 
2060; 3,846 cards 


Note 1. The Admiralty Computing Service is now at an end ; but the punched 
cards listed under [67] are available. 
Definitions: D? = + 6:7); D' = 6, — 
(If 4 is the interval between two consecutive arguments, then f(x» + ph) 
= fo = fot pD' + pD*.) 
Note 2. Fermi functions 22.68 


(2» + 2)!72 + 
A [ y! | (2p ro)?! ) T*(2s, + 1) 


= 


y=asW/p; W=VpP+1; = fa Al; 
a = 1/137.03 


z = atomic number of 8 emitter; z > 0 for electrons; z < 0 for positrons. 


A 
N 


; 
( 
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A = mass number of 8 emitter; = momentum of £ particle in mc units 
W = total energy in mc* units; ro = nuclear radius in h/mc units. 
Write to (4) for further details. 


Note 3. Tables relating to hydrodynamics 26.2 


w (x,y, = 1— (— (us) F 9,0) 
where 
F (n, y, 0) = (exp 4; + exp u2)/(1 + exp us) 
uy = — (m+ 4) (1 — = — + +9) 25 = 1) 5; 
us, = (n + 4) xx. See Bureau of Mines R.I. 4885 (July 1952) for further 
details. 


Note 4. The random digits were prepared by W. F. Brown, Jr., of Sun 
Physical Laboratory, Newton Square, Pa., with the aid of random digits of 
the RAND Corp. 
Corrections: [MTAC, v. 5, 1951] 

p. 190, Source 55(b), for sinh F — 1 read sinh F — F. 


p. 196, 7.1 Items 2 and 10 under this listing, (18)C : for x =0(.001)1.571; 
7D read x = 1(.001)1.228. 


p. 201, 10.0 Item 9, 1. 13. for 14S and 18S read 15S. 
p- 201, 11.6 Add the remark: “‘See also 28.3.” 


p. 203, 14.6 forx = 9(.1)10;y = 0(.1)10; 14D or 15D read x, y=0(.1)10; 
12D and 14D. (Extension previously in progress now com- 
pleted.) 


. 205, 17.1 Item 3, add ‘‘Available at (18).” 

. 209, 25.2 Tabulation of x? not available at (4). 

. 209, 25.3 for Cards 0(1)14,566 read Cards 0(1)13,000. 

. 209, 25.3 for 9.4247 read 1.5707. 

. 211, 28.2 Kepler’s equation: for E = M — esinE 
read M = E — esin E. 


Keypunching in progress. 
(a) At (24): I(x), 3, i, 3 


(b) At Point Mugu, under supervision of Donald Dufford, NAMTC, 
J,(x), for same values of » as those above. Thus both volumes of the 
NBS Tables of Bessel functions of fractional order are now being 
keypunched. 
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(c) At [24]: I(x), based on manuscript of BAAS now in print. 
(d) [24]: Spherical Bessel functions. Keypunched but not checked. 


(e) At NOTS, KennetH C. and CHarLes RICKER, China Lake 
Pilot Plant, China Lake, California. Bessel functions y,(x) and 
Y,(x), 2 < 20, BAAS, v. 10. 


Readers are again requested to review their keypunching loads, for 
possible keypunching of other basic mathematical tables during spare hours. 
Any one who can undertake a part of this work is requested to communicate 
with the undersigned. 

GERTRUDE BLANCH 


EVERETT C. YOWELL 
National Bureau of Standards 
Los Angeles 24, Calif. 


The Sieve Problem for All-Purpose Computers 


Introduction. The term all-purpose digital computer is often used to 
indicate a computing machine capable of performing the rational operations 
using addition and multiplication as basic functions of the arithmetic 
unit. Even when these operations are supplemented by some discrimination 
and “extract’’ commands, for dealing directly with the digits of numbers, 
there are a number of processes to which such a machine is not well suited. 
This includes even finite processes that are combinatorial in nature. Perhaps 
the most well-known of these processes is the sieve process. Although special 
equipment has been constructed for carrying out this algorithm! we shall 
not describe it here. Our purpose is to indicate how the all-purpose computer 
may be programmed to compete with sieve machines. 

Let us first state the general problem to be solved by the sieve process. 
Let m, mo, ---, m, be a set of k positive integers which, for the purposes we 
have in mind, may be assumed to be relatively prime in pairs. For each m; 
we consider ; distinct arithmetical progressions or linear forms in the 
variable x which we denote by 


We may assume that, for 7 fixed, the a;; are distinct non-negative integers 
less than m;. The problem is to find all integers N between given limits, say 


A<N<B, 


such that each N belongs to & arithmetical progressions. 

The number & is called the width of the problem, the k numbers m; will 
be called the moduli of the problem and the numbers a;; (j = 1 --+ ,) the 
admissible remainders for the modulus m;. The solution N, or even the 
number of solutions, is an exceedingly complicated function of the given 


parameters m;, a;;, A, B. There are two extreme cases which may be 
mentioned. 
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Eratosthenes Sieve Problem. In this case m, = 2, m, = 3, and in 
general m; is the i-th prime number, k is the number of primes not exceeding 
A = Bi, n; = m; — 1 and a, =j. This gives the famous sieve of Eratos- 
thenes and has for solutions (besides N = 1) all the prime numbers between 
Bt and B. In this sieve there is the maximum number of admissible re- 
mainders for each modulus. 

Chinese Sieve Problem. In this extreme case there is the minimum 
number of admissible remainders for each modulus, that is, 2; = 1 for all ¢. 
This problem is very old? and often “solved” by what is known as the 
Chinese Remainder Theorem. If A = 1 and B is the product of the moduli, 
then there is a single solution N which may be found by any one of a number 
of other nearly equivalent practical methods rather than by a bona fide 
sieve process. On the other hand, this type of sieve problem is often useful 
as a checking routine for a more general sieve setup. 

Quadratic Sieve. Midway between the two extreme examples we have 
examples in which each m; is approximately $m; so that, to put it roughly, 
any number N has, a priori, an even chance of belonging or not belonging 
to one of the arithmetical progressions of a fixed modulus m. The expected 
number of solutions N is therefore approximately (B — A)/2*. This kind of 
problem is the one most frequently met with and occurs in problems in- 
volving quadratic residues and congruences, binary quadratic forms, 
Diophantine equations of the second degree, etc.; hence the reason for 
calling this the quadratic sieve problem. 

Theoretical Aspects. It is theoretically possible to combine any two 
arithmetical progressions with relatively prime moduli 


Py = mx+ ay = mx + 
into a single one of the form 
P; = mmx + as 


where P; consists of all numbers common to P; and P2. Thus all numbers 
belonging to both the forms 


5x +4 12x +5 
comprise the arithmetical progression 
60x + 29. 


By proceeding in this way one may combine the & sets of arithmetical pro- 
gressions into a single set whose modulus is m = mym, --- m,. The number 
of arithmetical progressions in this one set is clearly = mym2 --- my. The 
ratio n/m is the density of solutions of the problem. The above process 
becomes intolerably unwieldy even for moderate values of k except when 
most of the m; are equal to unity. In the quadratic case, for example, a 
moderate value of k, like 20, which eliminates all numbers but one in a 
million, would lead to one set of at least 4-10" arithmetical progressions 
with a modulus of 4-10*". To produce and sort for size the 4-10" admissible 
remainders is clearly out of the question. In these cases it may be, indirectly, 
more practical to examine all integers N between A and B excluding each 
integer in turn for non-membership in one of the original sets of arithmetical 
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progressions. This is the procedure we speak of when we refer to a sieve 
process. 

Normalized Sieve Problem. In order to make a comparison of the 
effectiveness of different sifting processes we restrict ourselves to the case 
in which »; > 1 for all 4. That this is no real restriction is seen from the 
following consideration. If, for example, m,; = 1 so that we are looking for 
an N of the form mx + au, we may change variable from N to N, where 
N = mN, + au and thus eliminate P(x). This reduces the width of the 
problem by unity, introduces new constants a;; in the remaining k — 1 
sets of progressions, but leaves the other m; unaltered. The limits A and B 
are replaced by A/m, and B/m, so that the result is a fictitious speed-up 
of the sieve process by a factor of m. Proceeding in this way with any other 
cases of m; = 1 we finally come to grips with the real problem in which 
nm; > 1. To be sure, a case in which n; = 2 breaks up the problem into two 
parts each of which may be considered separately as a problem with an 
nm = 1. The result is two short problems instead of one long one. This dis- 
section of the problem is a favorite method with hand computers and 
accounts for the fact that fairly slow automatic sieves have considerable 
competition from hand methods. Of course the same dissection method is 
applicable to automatic sieves too, but often it is not practical because of the 
comparatively long set-up time involved. For these reasons, therefore, we 
consider as a normalized sieve problem, one in which n; > 2. To put the 
matter in another way, we consider that all integers N between A and B 
have, a priori, an equal chance of satisfying the conditions imposed by the 
problem. The effectiveness of a particular process will then depend only 
upon the speed with which it can dispose of the average candidate N and 
pass on to N + 1. 

Rates of Special Sieves. Without going into details it is desirable for 
purposes of comparison to mention existing sieves and to give an idea of 
their effectiveness. First of all we have hand methods. These involve the 
so-called movable strip method and its generalization, the stencil method. 
Here the range of natural numbers is represented by one or more sheets of 
paper ruled in squares, each square representing, by virtue of its position on 
the sheet, a definite integer. Strips of paper of length m; are punched with 
m; — n; holes, representing the inadmissible remainders modulo m;, and are 
moved carefully down the columns of the sheets and those square cells which 
are revealed by the punched holes are crossed out by the computer. After 
k of these strips have been applied, those cells which still survive are the 
answers to the sieve problem. Rates as high as 3 numbers per second are 
difficult to maintain over a long period of time. 

Electromechanical sieves have been built to canvass 50 numbers per 
second. These sieves can accommodate almost any reasonable sized moduli 
m; and handle problems of width k < 20. 

A photo-electric sieve has been built in 1932 which runs at 6000 numbers 
per second and an electronic sieve is under development which is designed 
for approximately 300000 numbers per second. These high-speed sieves are 
not equipped with high-speed output and are intended to be used on prob- 
lems whose progressions give restrictions comparable to or higher than those 
of quadratic type. It will be seen that such high speeds are not obtainable 
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with all-purpose computers but that the more modest speeds of the electro- 
mechanical type sieve can be achieved and even surpassed. 

Direct Attack by an All-Purpose Computer. Most all-purpose com- 
puters are equipped with a division command or at least a division program. 
A direct method of handling the sieve process would consist in dividing a 
trial value of N by each modulus m; in turn and then inspecting the re- 
mainder for membership in the set a;;. In the case of the Eratosthenes sieve, 
the inspection of the remainder is practically instantaneous as it requires 
only one or two addition times. If the number of moduli is large, for example 
large enough to make a list of primes above 10’, then the number of divisions 
in case N is a prime will be greater than 440. Speeds are such that 50 milli- 
seconds per division is not often realized. Allowing only 10 milliseconds per 
division, times of the order of several seconds are required to treat these 
cases. This does not compare favorably with even hand methods. However, 
in case the number of moduli is small a much more favorable speed prevails. 
A small Eratosthenes sieve is often incorporated into a problem in which a 
parameter » is supposed to have a prime value. In this way a large per- 
centage of composite p’s are eliminated, the remaining composites being 
taken care of later either from the output or by some other internal program. 
Beyond a certain point we get diminishing returns from each new modulus 
we introduce. 

In the case of a typical quadratic sieve problem the examination of the 
remainder would be a much longer routine but still, on the average, rather 
shorter than a division program. Sieves of width 25 would canvass numbers 
at the rate of three or four per second. This is still too slow to compete with 
much simpler equipment. Besides, the memory capacity of many machines 
would be exceeded in trying to accommodate the given remainders a,; in 
many typical problems. 

Binary Method of Representing a Sieve Problem. The fact that, as far 
as a given modulus m is concerned, a proposed number N is either rejected 
or accepted, makes it possible to treat the sieve with binary methods. The 
sieve problem may in fact be represented by a matrix of k rows and infinitely 
many columns. The infinite rows are periodic, the i-th row having a period 
of m;. The element situated in any column whose number is congruent to r 
modulo m; is 0 or 1 according as r is or is not one of the acceptable remainders 
a;;. The problem is then to find those columns which consist wholly of 
zeros. (Clearly the réles of 0 and 1 may be interchanged here.) The matrix 
will be referred to as the sieve matrix. 

For example, suppose that the problem is to find the least positive integer 
of the forms 


7x +2or3o0r4 
19x + 1 or 60r7 


fret lor4 


The matrix corresponding to this problem is 


t 

01101011010110--- 

10001111000111--- 
111001101111--- 
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The answer to the problem will be seen to be 16. In fact if we continue each 
periodic row two more steps beyond what is shown, we obtain a column of 
zeros. 

Machine Methods of Realizing the Sieve Process. Most of the high- 
speed digital computers use binary numbers and almost all of the decimal 
type machines use a binary code for decimal digits. The purely binary 
machines are the most natural ones with which to apply the following 
methods. The use of decimal type machines entails a considerable reduction 
in the width of the problem or else a certain amount of indirect use of the 
binary code for the decimal digits. In what follows we shall assume that we 
are dealing with a binary type machine, leaving the modifications necessary 
for the decimal machine to the ingenuity of the reader. The operations 
needed to carry out the sieve process clearly involve shifting operations 
and the extraction or discrimination of digits. This does not mean that 
addition, subtraction, and multiplication are not needed. In fact these latter 
operations are frequently involved in the execution of the former ones. 

Of the various possible methods to be used that one is to be preferred 
which gives the greatest speed in canvassing the columns of our matrix. 
The decision may be determined by the logical structure of the machine. In 
any case, there are two general steps to carry out: (a) the passage from one 
column of the matrix to the next and (b) the inspection of a column to see 
if it contains a non-zero element. 

First Method. The steps (a) and (b) can be carried out simultaneously 
by the simple process of doubling the binary number represented by each 
row. Since in many cases the period m; will exceed the number s of digits 
which the adder wili accommodate, it will be necessary to program this 
doubling as a ‘multiple precision’ process, taking account of overflow 
between the separate batches of s digits of the entire number. In particular, 
if overflow occurs in the first column there are two important consequences. 
First, this indicates that a 1 was standing in the first column and so the 
corresponding number JN is to be rejected. Secondly, this overflow must be 
sent to the other end of the period to perpetuate the periodic pattern of 
the digits. This operation has to be performed k times, once for each row. 
If at any time no overflow is obtained in the entire set of k operations, then 
we have a solution of the problem and the machine is instructed to print 
out the answer, continuing on or halting as desired. The entire process can 
be carried out using only addition with detection of overflow, and involves a 
certain amount of “‘tallying’’ and modification of commands. The method 
has the disadvantage of requiring k multiple precision operations in passing 
from N to N + 1. When £ is large (20 or 30 in some typical cases) and the 
m’s large, the major portion of the time is spent in keeping the periods up- 
to-date. However, for small & and m, the method is both simple and fast. 
In fact, its use is recommended in problems in which no actual sifting is 
involved. Suppose, for example, that we have a problem involving one or 
more parameters and these are given certain irregularly spaced values such 
as for instance 


x = 1, 3, 6-9, 12, 15, 18-22, 25. 


Instead of storing all these numbers in as many memory positions, this 
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information may be stored in only the one word 
1010011110010010011111001, 


in which the r-th digit from the left is 1 or 0 according as r belongs to our 
set or not. By successive doubling and inspection for overflow the machine 
can inform itself as to which values of x it should consider. 

Second Method. In order to eliminate the lengthy multiple precision 
processes of the first method one may proceed as follows. The periodic 
digits of each row of our matrix are stored as before in ‘“‘words”’ of s digits. 
The minor of & rows and the first s columns of the original matrix are now 
duplicated in a special part of the memory where the operation of doubling 
and testing for overflow is performed as in the previous method. However, 
when overflow occurs the fact is recognized by the machine but the digit 
is not kept track of as before. After s doublings the information in the 
original k by s matrix is all used up and the matrix now consists wholly of 
zeros. 

The next process consists of circulating the original matrix by moving 
the elements of each row s places to the left. Since s binary digits is a word, 
this process is simple and fast. Only one complication arises, that of dis- 
posing of the first word in each row. If for a particular modulus m; we have 


then the first s — r, digits of this first word w; belong in the penultimate 
word of the shifted matrix while the last r; digits form the last word. This 
disposal of w; can be accomplished by use of a product command. The 
machine computes the exact two word product of w; by 2*-". The “most 
significant” word of this product is thus added to the penultimate word of 
the shifted matrix and the least significant word is placed at the end. 

One disadvantage of the second method is the fact that all k rows of 
the matrix have to be treated, even though a rejection of the value of NV 
may result from one of the first few rows. This early rejection of N is a 
highly probable event in the case of a quadratic sieve problem. In fact the 
probability of getting a rejection after h trials is 1 — 2-*, so that only one 
N in a thousand is apt to pass the test imposed by the first ten rows. More- 
over, in many problems some moduli are much more restrictive than others 
and if the rows of our matrix which correspond to these moduli are put first, 
the probability of early rejection is very considerably increased. In order 
to exploit this possibility we propose another method. 

Third Method. This method differs from the previous one in its more 
rapid treatment of the minor of & rows and s columns. Instead of using the 
detection of overflow we employ the more elaborate “‘extract’’ command. 
This command, which varies slightly from one binary computer to the next, 
enables the machine to extract from a given binary number those digits 
which occupy specified positions. These positions are specified by a number 
called the extractor. As extract is performed in the SWAC, for example, 
those digits of a given number, called the “‘extractee,” are selected which 
correspond to the zero digits of the extractor; all other digits are made zero. 
Thus if the extractee is 


11101010001--- 
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and the extractor is 
20408 


the extracted number is 
10000010000---. 


In other words, the extracted number has for its i-th digit the product of 
the i-th digit of the extractee by the complement of the #-th digit of the 
extractor. 

Let us now consider the minor of & rows and the first s columns of our 
sieve matrix. Starting with the extractee consisting of s 1’s and using the 
first row of our minor as extractor, the extracted number records by its 
1’s those columns whose first row elements are zero. Using this extracted 
number as a new extractee and the second row as extractor, the extracted 
number now records those columns whose first and second elements are 
both zero. Continuing this process and comparing the successive extracted 
numbers with zero (that is the number consisting of s digits 0) after each 
extraction, we very soon (at least in the typical case) arrive at an extracted 
number which is zero. At this point the entire minor is thrown away and a 
new minor is created by circulating the rows of our matrix as described in 
the second method. If, however, our minor contains one or more columns 
of zeros, a non-zero number will result from all k extraction operations with 
digits 1 in the positions corresponding to the solution or solutions of the 
problem. To find these positions we may proceed successively to double the 
corresponding number and test for overflow as in Method 1. An obvious 
method of tallying will produce the solution or solutions which the extraction 
process has detected. 

Comparison of Methods. On the basis of speed, a comparison of the 
three methods for the case of the quadratic type sieve problem may be made 
as follows. 

Let 


T = the average time required to dispose of a single number NV and 
pass on to NV + 1, 

s = the number of binary digits in a machine word, 

h = the average number of words for each modulus of the problem, 

k = the number of moduli, or width, of the problem, 

a = the addition time of the computer, that is, the time required to 
receive and add two numbers and to dispose of their sum. 


For the three methods rough approximations to T may be given as 
follows: 


First Method: T = Shka 
Second Method: T = 5k(1 + 2hs“)a 
Third Method: T = [10kh + 6(1 + loges) ]s—ta 


As an example of what may be done on the SWAC for which s = 36 
we take a typical case of h = 2 and k = 17. The addition time for the 
SWAC being 64 microseconds, the three values of T turn out to be 10880, 
6044, and 672 microseconds respectively. 


Nn 
F 


re 
ti 
1 
7 
1! 
n 
‘ 
‘ 
a> ‘ 
” 


n, 


to 


1e 
0, 


THE SIEVE PROBLEM FOR ALL-PURPOSE COMPUTERS 13 


The above formulas for T do not take into account certain occasional 
routine operations such as resetting tallies to zero, the restoration of modified 
commands, the decimal conversion and printing of the solutions, etc. To 
take these activities into account may lengthen the time by as much as 
10 percent. In actual tests on the SWAC, for instance, the third method gave 
T = 695 microseconds. It should be pointed out that 4 must be at least 2 
in the second and third methods since a complete minor of s columns is 
needed. 

Operating Suggestions. We conclude with a few remarks for the benefit 
of the programmer and operator in dealing with sieve problems. In the 
normal case the expected output is only a few numbers N occurring in un- 
predictable places between A and B. It is evident that the machine must be 
in absolutely perfect operation during the run in order to be sure that no 
solution has been overlooked. Absolute precision in the circulation of the 
rows of the sieve matrix is, of course, of paramount importance. The ac- 
curacy of this operation may be checked at the end of the run by printing 
out the entire matrix and inspecting the result. In order to check the machine 
during operation the programmer may deliberately insert one or more 
“traps” or precomputed extraneous solutions. It is not difficult to find by 
hand methods a number N>» which satisfies most of the requirements imposed 
by our matrix. To make N> satisfy all the conditions we may deliberately 
weaken the remaining ones by changing a 1 to a 0 in each of the outstanding 
rows. This weakening, especially when applied to large moduli, will not 
change perceptibly the overall exclusion ratio but may introduce one or 
more unexpected extraneous solutions which, however, are easily recognized 
when they occur. This is a very cheap price to pay for the confidence that 
one gets from seeing the machine deliver the one or more predicted solutions 
No on schedule. 

In case the running of the problem is interrupted either by machine 
failure or by another problem of higher priority, there are three procedures 
available. 

First, one may start the problem over from the beginning. This is 
especially advisable if the interruption occurs early in the run. Secondly, 
one may make an obvious translation in the variable VN and compute by hand 
a new sieve matrix. Of course, this has to be done with considerable care but 
in case the machine is being used on another problem, there may be time 
to do this. In running very long problems it may be wise to prepare in 
advance several matrices from which a new start may be made in case of 
interruption. This procedure, however, is almost as much work as breaking 
up the problem into parts as described above, and this latter procedure 
always pays off in the consequent speed-up of the problem. 

A third procedure is available in case there is still some memory space 
unused. This method is simply to make the machine compute its new starting 
matrix itself. This program may be based on the First Method, using 
detection of overflow to circulate each row of the original matrix to its new 
position. Of course the amount of this circulation depends on the new 
proposed starting point and is different for each modulus. If, instead of 
starting with N = 1, we wish to recommence with N = N,, the row cor- 
responding to the modulus m must be circulated rj]—}1’steps where’ r is the 
remainder on division of N, by m. The program for carrying’ this out is set 
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up with N, as a free parameter which can be typed into the machine as 
occasion demands, no further information being needed. This elaboration 
of the program makes it possible to operate a lengthy sieve problem as a 
backlog workload without the need of a specially trained operator. 
D. 1. 
1D. H. Lesamer, “The mechanical combination of linear forms,’ Amer. Math. 
Monthly, v. 35, 1928, p. 114-121, “‘A photo-electric number sieve,’’ Amer. Math. Monthly, 


v. 40, 1933, p. 401-406, ‘‘A machine for combining sets of linear congruences,” Math. 
Annalen, v. 109, 1934, p. 661-667. 


pin in = Dickson, History of the Theory of Numbers, v. 2, Washington, 1923; New York, 
» 


The Use of Large Intervals in Finite- 
Difference Equations 


In a recent article Stk RICHARD SOUTHWELL’ has challenged the general 
theory of finite differences and in particular the use of it in connection with 
the solution of differential equations by relaxation methods.” Opinions differ 
on the method of treatment of a numerical differentiation formula such as 


1 


where h is the interval between pivotal points and where 6?"yo is the 2n*® 
central difference of ‘yo. This equation is replaced for the purpose of numerical 
solution by its equivalent 


= (yi + y-1 — 290) + A (yo), 


where 41, Yo, and y_; refer to the values of y at xo + h, xo, and x» —h 
respectively and 


1 1 


I advocate the use of as large an interval as possible (consistent with con- 
vergence of the finite-difference equations and convenience of computation) 
calculating A(y) and incorporating it into the computation. Southwell 
prefers to use such a small interval that A is negligible. He states: ‘Accuracy 
is not predictable from quantities computed from a finite-difference approxi- 
mation unless the interval is small enough to justify belief in the convergence of 
the Taylor series: in general the radius of convergence, though it exceeds one 
or two of the smallest intervals that are practicable, will not exceed many 
such intervals, and consequently what have been claimed as closer approxi- 
mations may in fact be less accurate.” 

The purpose of this note is to point out that, in the writer’s opinion, this 
statement is unduly cautious and that its conclusion cannot be true if finite- 
difference equations are properly used. 

In the first place it is not clear what role the Taylor series plays in finite- 
difference theory. If the Taylor series is not convergent, care must certainly 
be exercised ; but finite-difference methods do not always break down in this 
case.* The point to be emphasized, however, is that an examination of the 
differences will tell us, in all practical cases, whether the finite-difference 
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equations we wish to use are valid and how many differences are significant 
at the particular interval used. It is stressed by Fox? that in all the finite- 
difference formulae used the differences must be convergent and that other- 
wise a smaller interval must be used. 

By ‘“‘convergence of differences’’ I mean that the differences of some order 
should oscillate about zero, with a maximum amplitude determined in the 
usual way by the binomial coefficients of that order.* The tabulated function 
can then be regarded as a polynomial, slightly perturbed by rounding errors 
of not more than half a unit at any point. The degree of the polynomial is the 
same as the order of the last significant difference, and the function can be 
interpolated, integrated, and differentiated by finite-difference formulae 
based on this polynomial representation, the error of the process being of the 
order of the first neglected term in the finite-difference equation. 

The degree of the approximating polynomial varies with the number of 
figures to which the given tabular entries are rounded off and with the magni- 
tude of the interval between successive pivotal points. The function 
(1 + x?)-, for example, can be represented near the origin, at an interval of 
.1, by a quadratic to two decimals, a sextic to four, and so on. In each case 
interpolation and integration can be performed effectively to the same pre- 
cision as that of the tabular entries, though derivatives suffer from the loss 
of significant figures in successive differences. In particular no estimate can 
be given for derivatives of the true function of orders greater than the 
degree of the approximating polynomial. 

The fact that the true higher derivatives may be large and ever increasing 
as for example in the function (1 + x*)-, is of no significance with regard to 
the use of finite-difference formulae as applied to the approximating poly- 
nomial. If the derivatives increase fast enough such formulae may be asymp- 
totic, but again the error is of the order of the first term neglected. This is 
analogous to the treatment, in a recent paper by Fox & MILLER? , of a func- 
tion whose Taylor series had zero radius of convergence at the origin; the 
series was in fact asymptotic. 

Southwell does not use differences and bases his argument on the fact 
that increased accuracy is not guaranteed by increasing the order of the 
approximating polynomial, citing as one example the interpolation at unit 
interval of the function (1 + x*)~'. The first table gives the differences of this 
function, tabulated to three decimals at unit interval (the central differences 
at x = 0 are filled in from consideration of symmetry). Examination shows 
that the differences are diverging rapidly, that interpolation formulae can 
hardly give even one-figure accuracy in the function near x = 0, and that a 
much smaller interval is necessary. 


x 10% x 10% 

0 1000 — 1000 2400 12 6897 285 
—500 1200 —1015 -16 

1 500 200 —1200 13 5882 209 22 
—300 0 — 806 —54 

2 200 200 —141 14 5076 155 16 
—100 —141 —651 —38 

3 100 59 15 4425 117 12 

—41 —534 —26 
4 59 16 3891 91 
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The divergence of differences near x = 0, however, does not prevent the 
accurate use of finite-difference formulae at points sufficiently remote from 
this point. Central-difference formulae, for example, are based on polyno- 
mials whose origin is at the particular pivotal point considered, and the func- 
tion (1 + x?) can be quite correctly interpolated, differentiated, or inte- 
grated to at least four figures, at unit interval, in the range shown in the 
second of the tables. 

Turning now to the rest of the quotation, the second part suggests that 
the largest interval at which the full finite-difference equation can give 
accuracy will never greatly exceed an interval for which its first term alone is 
adequate. This will depend, of course, on the particular function involved, but, 
even if the ratio of intervals is only two, the use of A is worth-while not 
only for economy of effort but also from the point of view of accuracy. 

Computers know that the use of small intervals in step-by-step integra- 
tion is dangerous, the tendency for rounding-off errors to accumulate being 
roughly proportional to the square root of the number of intervals used. In 
relaxation methods, similarly, a small interval involves a large number of 
simultaneous equations, and their accurate solution is a matter of consider- 
able difficulty ; often one or even two extra figures have to be retained in the 
coefficients and residuals. Practical views on this have been reinforced by a 
recent theoretical paper by Topp.°® 

In illustration of these points let us try to evaluate to four decimals the 
pathological function (1 + x*)—! from the facts that it satisfies the differential 
equation 


(1 + y’ +2 (14+ y = 0, 


is symmetric about x = 0, and has the value y ($) = .8000. The finite- 
difference equation is given by 


(1 + 2x0 fo) + y-1 (1 — 2x0 fo) — yo (2 — 2hfo) + A(yo) = 0, 


where 


and 
f=h(A t+ x’) 


At an interval h = .25 there are only five pivotal points, giving insuffi- 
cient information about the differences. The next smallest convenient inter- 
val is .1 (assuming we wish to have x = 0 as a pivotal point), and at this 
interval we find the following first approximation y, with A neglected. 
Extra figures have been retained in the residuals and coefficients to ensure 
that this approximation is everywhere an accurate solution (to four figures) 
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to the given simultaneous equations. The evaluation is as follows: 


Error in 

x A R 104y® 

0 10020 — 200 22 —1.8 10001 20 
—100 11 

9920 —189 21 —-1.9 —2.1 9902 19 
—289 32 

2 9631 —157 12 —1.5 —1.5 9616 16 
—446 44 

a 9185 —113 2 —1.0 — 8 9175 11 
—559 46 

A 8626 —67 (-—9) — 2 — 2 8621 5 
—626 (37) 

$ 8000 8000 0 


The differences look quite satisfactory, fifth and higher orders being 
apparently negligible, and the recorded values of A include contributions 
from third and fourth differences only. If these values are now inserted in 
the finite-difference equations, the remaining residuals are as shown in the 
column headed R, and their liquidation gives a better approximation y™, 
whose maximum error is one unit. The skilled computer would notice further 
that the fifth and sixth differences, though small, have a definite flow and 
are not yet oscillating about zero. Their total contribution is only about —.2, 
but their effect is to reduce the maximum error in y™ to half a unit in the 
last figure. 

This work involved the solution of nine simultaneous equations with 
two sets of residuals and with a little intermediate computation. If A were 
ignored and the computer solved 19 linear equations at an interval of .05, he 
would still have maximum errors of five units. 

It is very doubtful whether intuition alone could here decide the question 
of the size of the “ultimate interval” for which A is negligible, and, if the 
computer, still ignoring A, wished to make sure of his result by taking an 
interval of .025, he would still have errors of one unit, having solved no less 
than 39 linear equations in addition to computing the coefficients thereof. 
The labour of this process is excessive; the accurate solution of these 39 
equations would involve keeping at least two extra figures in the coefficients 
and residuals. 

I therefore maintain, with regard to the end of the statement, that a 
computer familiar with both the theory and practice of finite differences will 
never be in danger of claiming more accuracy than he has achieved. Either 
he has got the desired accuracy, or he will know that he has not got it (and 
perhaps cannot get it, as for example near a discontinuity or singular point). 
No question of prior knowledge of the convergence of the Taylor series 
arises, and such knowledge has certainly not been sought in any of my 
published examples. Examination of the differences, after a tentative solution 
has been found, provides the necessary information and is of course essential. 
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In this paper I have concentrated on functions of one variable, but the 
remarks apply with greater force to partial differential equations, particu- 
larly in the desirability, both for the sake of accuracy and a minimum of 
labour, of using conveniently large intervals. 


L. Fox 
National Physical Laboratory 


Teddington, England 
and NBSCL 

1 Numerical Methods of Analysis in Engineering. Edited by L. E. GRinTER, Macmillan 
Compete, New York, 1949, chapter 4. 


Fox, “Some improvements in the use of relaxation methods for the solution of 
ores. and partial differential equations,” Royal Soc., London, Proc., v. 190A, 1947, 


3L. Fox & J: C. P. Miter, sees for large arguments. The exponential inte- 
gral,” MTAC, v. 4, 1951, . 163-167. 


for example, MILLER, “Checking by differences—I,” MTAC, v. 4, 1950, 


p. 3-11 
‘7. T opp, “The conditions of a certain matrix,’’ Cambridge Phil. Soc., Proc., v. 46, 
1950, p. 116-118. 


Monte Carlo Matrix Inversion and 
Recurrent Events 


1. Introduction. Recently Wasow! has given a necessary and sufficient 
condition that one of two unbiased estimators of the inverse element of a 
given matrix has a smaller variance. Using the theory of recurrent events, 
we extend a result by FELLER? in order to generalize and reinterpret Wasow’s 
condition. To do this let us consider a simple discrete Markov process with 
a finite number of states. Of these m + 1 states denoted by 0, 1, 2,---,m, we 
prescribe that the state named ‘‘0”’ is the only death state or sink—in the 
sense that the random walk ends when this state is reached. Let p;, be the 
one-step transition probabilities for i, k = 1, ---, m, and let 


po = 1—-— pix for «= 1,---,m. 
k=1 


Further, we assume that each pio > 0. When x > 0 define p{? as the proba- 
bility that in a random walk starting at state 7, state k is visited on the 
n-th step; and define p{? = 8:x, where 5,4 is the Kronecker delta. 

Now if (qx) = (din — pix) is the m by m matrix whose inverse (q**) is 
desired, we may estimate each element g** by Monte Carlo methods as in 
Wasow’, ForsyTHE & LEIBLER’, and Curtiss‘, since 


qi* ae PR. 


Two estimators have been investigated. The first estimator is the sample 
mean of the random variable s;, where six = pio! if the random walk having 
started at state 7 stops just after visiting state k, and s;, = 0 if otherwise. 
The second estimator is the sample mean of the random variable v;, where 
vx is the number of visits to state k in a random walk having started at 
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state i. Now it is known that both these estimators are unbiased, that is, 


(1) E(six) = 
and 
(2) E(vix) = q* 


where E is the expected value operator. Further, the second moments of 
these random variables can be expressed directly in terms of the g** and 
Peo. We have 


(3) E(siz) = 
and 
(4) E(viz) = (2g** — 


Relations (1) and (3) were proved by ForsyTHE & LEIBLER’, and relations 
(2) and (4) were derived independently by Curtiss‘ and by the author. 
Now, Wasow! showed that 


Nix 


(3’) E(siz) = fori #k 
and 
(4’) E(vn) = + Aas) fori #k 


(1 — Age)? 


where he defined X, as ‘‘the probability of going from state i to state k 
without passing through state k on the way; i.e., Aix is the total probability 
associated with all paths connecting state i and state k, all intermediate 
points being different from state k.”” Thus, we have before us two different 
interpretations of the second moments of s;, and 0, in the case 7 ¥ k. 

2. Recurrent Events. Using the concept of recurrent event as developed 
by Feller? we now extend a result of Feller? to establish a general theorem 
which shows the fundamental probabilistic relation between the two inter- 
pretations for the moments. When n > 0 define rf? as the probability that 
in a random walk starting at state i, state & is visited for the first time on the 
n-th step, and define rf? = 0. Let 


n=0 
so that 7; equals the probability that in a random walk starting at state 7 
state k is eventually visited. 
THEOREM: — = rig*. 


Proor: By considering the mutually exclusive and exhaustive cases it 
follows that 


= ---+ Ao? forn>0 
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and since r{? = 0 and p§? = 5% that 
(*) = forn>o. 


j=0 
Next we define 


Pix(u) = 
as the probability generating function for the sequence {p%?}, and define 
Rix(u) = 


as the probability generating function for the sequence {r$?}. Hence by a 
basic property of generating functions we have by (*) 
Pix(u) — = Ru(u)Pix(u) for all u. 


Putting uw = 1 we get 


However, this is 


q* — = rag* 
which was to be proved. 
As corollaries we get the explicit relations 


qi* ome Siz 
an 
(6) q 


since by definition 

g** > =1>0 
and by assumption 

1 — rex > Peo > O. 


Thus we have a direct derivation of the general relation between g‘* and 
riz for all and k. 

As a matter of interest, while working with Wasow, the author first 
derived relations (5) and (6) indirectly by equating expression (3) to (3’), 
and expression (4) to (4’). Thus the relation of g* to \;, was found. Next it 
was seen that ,x is, in fact, just 7;z. An informal verification is as follows. 

The probability that a random walk initially traverses a given path is 
equal to the probability that the walk takes that path times the probability 
that the walk thereafter follows some path from the class of possible paths 
remaining. But the latter probability is one. Therefore \;, may be replaced 
by Vike 

We may now write the second moments of our two random variables in 
terms of riz. 


(7) B(sh) = 
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and 


(8) E(viz) = 


+ rie 
— 1 — rer 


Thus, a comparison of (3) and (4) ue 


(9) E(six) < E(viz) if and only if < 2g** —1 
and similarly, (7) and (8) give 
1 1 1 + ree Tkk 
(10) < E(viz) if and only if — < ——. 
Pro 1 — re 


This, of course, is Wasow’s result expressed in terms of r4, rather than Axe. 
Since 74, is usually unknown, but riz > Pix; (10) implies 


(11) Else) < if 
These recurrent events associated with this stochastic process have 
permitted the derivation of further results which will be treated in a subse- 
quent paper. 


H. P. E>DMUNDSON 

University of California 
at Los Angeles 
198 W. ah ena “A note on the inversion of matrices by random walks,"” MTAC, v. 6, 

p. 
¥ 50 An Introduction to Probability Theory and its Applications, v. 1, New 

or! 

3G. E. Forsyrue & R. A. LerBLER, “Matrix inversion by a Monte Carlo method,” 

MTAC, v. 4, 1950, p. 127-129. 
4J. H. Curtiss, “Sampling methods a: eae to —— and difference equations,” 

IBM Seminar on Scientific Computations, 
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1042[A].—C. E. Fr6BerG, Hexadecimal Conversion Tables. Lund, 1952, 20 
p., 21.7 X 14.6 cm. 


These tables are designed to assist in the coding of problems for binary 
computers. There are four tables. Table 1 is a table of the integers 
1(1)2!9(2*)2" in both decimal and hexadecimal notation. Table 2 gives the 
hexadecimal equivalents of ”-10-**, m = 1(1)100, k = 1(1)8, correct to 13 
hexadecimals. Table 3 gives the hexadecimal equivalent of 50 frequently used 
constants correct to 16 hexadecimals. Finally Table 4 gives the decimal 
equivalents of numbers of the form 


n-16-* m = 1(1)15,k = 1(1)13. 


Results are given to 16D. The letters A B C D E F are used as hexadecimal 
notation for the numbers ten through fifteen. A few examples of the use of 
the tables are given in the introduction. This handy booklet should find its 
way into many a coding room. 
D. H. L. 


| 
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1043[F].—M. Kraitcuik, Introduction a la Théorie des Nombres. Paris, 
1952, viii + 202 p., 16.5 X 25.4 cm. 


The material for this work is taken largely from the author’s Théorie des 
Nombres,! v. 1 and 2, and Recherches,? v. 1. There are numerous small tables 
which may be reported as follows: 


(1) p. 


(8) p. 


(9) p. 


(10) p. 


(11) p. 


2. Factorization of + 1 + pipe---pn, where p, is the m-th prime, 
for n < 16. , 


. 4. Factorization of 2” + 1 as far as known in 1947. 
. 5. Factorization of + 1 + m! for n < 22. 
. 12, 13. Complete factorization of 2" —1 for m odd and 


ni = 1(2)99, 105, 107, 111(2)117, 123, 127, 135. 


. 38, 39. Factorization of 2"-+ 1 for odd m < 135 and for even 


n < 148 with some gaps. 


. 40, 41. Factorization of 10* + 1 for m = 1(2)31 and m = 2(2)20, 


24, 30, 36, 50. 


. 55-62. Tables of power residues. These give the least positive or 


absolutely least residues of k-th powers for prime moduli < LZ 
for the following values of k and L 

& & = & F 

L = 200, 200, 200, 422, 242, 632, 1000 

124-129. Tables of linear divisors of x? — Dy*. These tables are 
in three parts. The first part gives the arithmetical progressions 
in which p must lie if D is to be a quadratic residue of p for all 
square free numbers D less than 35 in absolute value. The 
second part gives this information for0 < D = 4k + 1 < 241. 
The third part is for 0 > D = 4k + 3 > — 239. These con- 
tain two errata, carried over from the corresponding tables in 
Théorie des Nombres, v. 1: 

p.125 D=157 for 107 read 109 
p.126 D=193 for 155 read 129. 

165-177. Tables of x (mod p) in x? — y® = N. These are for all 
primes p < 67 and are taken from Théorie des Nombres, v. 2, 
p. 156-166. 

182. Table of a cyclotomic quadratic form. This table gives the 
coefficients of the polynomials X, Y, such that 

. | On(x) = 4k +1 

— mxY,2(x) = n= 4k+3 
for all square free m < 35 and also 39, 42, and 51. It is taken 
from Recherches,’ v. 1, p. 88, with two errata. The coefficients 
of X29 should be 1, 15, 33, 13, 15, ---, 15, 13, 33, 15, 1. 

184-186. Factorization of 2‘"** + 1. The table extends to all 
4n + 2 < 162 as well as 170, 174, 182, 186, 190, 198, 210, 222, 
234, 250, 258, 270, and 330, with some results partially incom- 
plete. For m = 94, 114, and 150, results are put in the wrong 
columns of the table. 

D. H. L. 


1M. Kraitcuik, Théorie des Nombres. v. 1, Paris, 1922, ix + 229 p., v. 2, Paris, 1926. 
2M. Krartcuik, Recherches sur la Théorie des Nombres. v. 1, Paris, 1924 
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1044[F ].—A. Natucct, “Osservazioni sul problema di Fermat,” Un. Mat. 
Ital., Boll., s. 3, v. 6, 1951, p. 245-248 
This paper contains a short table of 


F,(x, u,v) = (x + u+ ov)" — {x™ + (x + u)*} 


forv = 1,u = 1,2andn = 3, 4, 5. The variable x ranges over integers until 
F becomes negative. Two or three negative values are given and x does not 
exceed 13. 

D. H. L. 


1045[F ].—Giusepre PatamA, “Numeri primi e composti contenuti nella 
forma 1848x? + yy? dell’ intervallo 11000000 — 11100000,” Unione Mat. 
Ital., Boll., s. 3, v. 7, p. 168-171, 1952. 

The author uses a recent list of primes of the eleventh million’ to find all 
the primes of the form 1848x? + y*® between the above limits together with 
values x and y. There are 202 primes. This is similar to a list above 10 million 
given by CUNNINGHAM & CULLEN.? The author also lists 121 composite 


numbers of this form. 
D. H. L. 
1J. P. Kuk, L. Potetti & R. J. Porter, Liste des nombres premiers du onziéme million. 
Amsterdam, 1951 [MTAC, v. 6, p. 81, errata p. 163 


1901.53 or CUNNINGHAM & J. CULLEN, idoneal primes,” Br. Asso. Adv. Sci., Report, 
Pp. 


1046[H ].—P. T. SMoLiAxov & A. N. Hovanskil, “K resheniu algebraiches- 
kikh uravnenii 3-1 stepeni’’ [On the solution of algebraic equations of the 
third degree], Kazan Filial. Akad. Nauk. SSSR., Jzvestita Fiz-Mat. 
Tekhn. Nauk, v. 1, 1948, p. 85-92. 

The authors take as a reduced form the cubic 


y+ Py+1=0. 


On p. 87 is a 2D table giving one real root y corresponding to 271 given 
values of P between — 10 and 10. 
D. H. L. 


1047[K ].—G. E. ALBert & R. B. Jounson, “‘On the estimation of central 
intervals which contain assigned proportions of a normal univariate 
population,”’ Annals Math. Stat., v. 22, 1951, p. 596-599. 


Let Y be a random variable with the cumulative probability function 
F(y) = (0? f exp[— (u — m)*/20* ]du. The statistic considered is 
A(¥,s;d) = + ds) — F(P — ds) 


N N 4 
where = Y,/N,s = (Y; — ¥)2/(N 


Wivks has shown! that if for given p, 0 < p < 1, one sets 
Apw = tel (N + 1)/N}}, 
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where f, is Student’s ¢, then the expectation of this statistic is 
=1-—p?. 


The problem treated in the present paper is: for given p, a, di, d2, such that 
0<p<i, 0O<a<1, O<d, =1-—p; 0<d,=p), find the smallest 
sample size N for which one has 


Prob. {1 — p — d; S A(Y,s;\pw) S1— p+ ds} Za. 


A table of such smallest sample values is given for p = .01, a = .95, .99, and 

for p = .05, .25, .50, a = .80, .95, .99, in each case for (di, dz) = (.075, .05), 

(.05, .05), (.025, .025), (.035, .015), (.05, .01), (.025, .01), (.02, .01), (.01, .01). 
Z. W. BIRNBAUM 

Univ. of Washington 

Seattle, Wash. 


1S. S. Wirxs, “On the a of sample sizes for setting tolerance limits,’ 
Annals Math. Stat., v. 12, 1941, p. 91-96. 


1048[K ].—D. H. Buarte, “‘A note on the significance levels of the distribu- 
tion of the mean of a rectangular population,”’ Calcutta Stat. Assn., 
Bulletin, v. 3, 1951, p. 172-173. 


For means of samples of » drawn from a continuous rectangular universe 
on the interval (0, 1) the author gives the 5%, 2.5% and .5% points to 4D for 
nm = 1(1)16. He remarks that for m = 16 the corresponding percentage points 
for the normal distribution function approximation are correct to 4D. 

©. 


1049(K ].—J. M. Cameron, “Tables for constructing and for computing 
the operating characteristics of single-sampling plans,’ Industrial 
Quality Control, v. 9, 1952, no. 1, p. 37-39. 


Let the probability of acceptance of a lot be denoted by P(A) and the lot 
fraction defective by ». For a single-sampling plan with sample size m and 
allowable number of defectives c, Table 1 provides 13 paired values of » and 
P(A) which can be used to plot the respective operating characteristic (OC) 
curves. The entries in Table 1 are values of mp which are given to 3D for 
c = 0(1)49 and for P(A) = .995, .990, .975, .950, .900, .750, .500, .250, .100, 
.050, .025, .010, .005. 

Table 2 is designed for constructing single-sampling plans whose OC 
curve passes through two points (f:, 1 — a@) and (2, 8) where ; is the frac- 
tion defective for which the risk of rejection is to be a, and where pz is the 
fraction defective for which the risk of acceptance is to be 8. The entries in 
Table 2 are values, given to 3D, of the ratio p2/p1 which are associated with 
three sets of paired values of a and 8; namely, (a = .05, 8 = .10), (a = .05, 
B = .05), and (a = .05, 8 = .01). Corresponding to the selected value of 
p2/p1 are values of mp; and c. The sample size is determined by dividing 
np; by p: and the acceptance number is read directly from the table. 

Since the tables were computed from the Poisson series, they are exact 
only when the probability distribution of the number of defectives in a 
sample of follows the Poisson law. However, for most practical cases, the 


tab 
Uni 
Mir 
105 
to 
Ch 
rat 
99 
Sti 
W 
x 
9 
I 
§ 
’ 


RECENT MATHEMATICAL TABLES 25 


tables will give satisfactory approximations when the distribution of de- 
fectives is binomial. 

G. W. McELRatTH 
University of Minnesota 
Minneapolis, Minnesota 


1050[K ].—J. H. Cuunc & D. B. DeLury. Confidence Limits for the Hyper- 
geometric Distribution. Univ. of Toronto Press, Toronto 1950, 72 leaves, 
22.5 X 30.9 cm. Price $2.25. 


The charts in this book for the hypergeometric distribution are similar 
to the ones for the binomial distribution given by CLoppperR & PEARSON.' 
Charts are given for population sizes 500, 2500, and 10,000, for sampling 
rates 5%, 10% (10%) 90%, and for confidence coefficients 90%, 95%, and 
99%. The methods of constructing the charts are given explicitly, as well as 
means for performing interpolation and extrapolation to cases not included. 

I. R. SAVAGE 
National Bureau of Standards 
Statistical Engineering Laboratory 
Washington, D. C. 


1C. J. CLoprer & E. S. Pearson, ‘‘The use of confidence and fiducial limits applied to 
the case of the binomial,’ Biometrika, v. 26, 1934, p. 404-413. 


1051[K ].—A. Hap & S. A. SinKBAEK, “‘A table of percentage points of the 
x*-distribution,”’ Skandinavisk Aktuarietidskrift, v. 33, 1950, p. 168-175. 


The present table gives values of xo? to 3D for which P(x? = x0") in which 
x? has f degrees of freedom for P = .0005, .001, .005, .01, .025, .05, .1(.1).9, 
.95, .975, .99, .995, .999, .9995 and f = 1(1)100. 


1052[K ].—H. Levene, “‘On the power function of tests of randomness based 
on runs up and down,” Annals Math. Stat., v. 23, 1952, p. 34-56. 


In the set of observed values, x1, x2, --*, x, of a continuous statistical 
variable X, let B* be the sequence of signs (+ or —) of the differences (xj: 
—x,) fori = 1, +--+, — 1. A sequence of successive + (—) signs not imme- 
diately preceded or followed by a + (—) sign is called a run up (down). Let 


1 
s be the number of runs up and let E’(s) = lim . E(s), o(s) = lim : o*(s), 


where E(s) is the expected value of s and o°(s) is the variance of s. Also let k 
be the total number of + signs in B*, with E’(k) and o’*(k) defined analo- 
gously to E’(s) and o’*(s). s and k are members of a class of u-run statistics 
which may be used for tests of randomness of the sample values, other mem- 
bers of which are also discussed by the author. 

One of the alternatives to randomness, precisely defined in this paper, 
is that x; obeys a normal distribution law with mean i# and unit variance, 
i = 1, ---, m. For this case E’(s), 1 — E’(k), o(k), and o’(k) are given to 
6D for the first three and to 3D for the fourth for 8/v2 = 0(.1)3(.2)4 in 
Table 1. An interesting graph of E’(s) for the linear trend considered com- 
pares the results for normal and rectangular universes. 

L. A. AROIAN 
Hughes Aircraft Company 
Culver City, California 
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1053[K ].—F. H. C. Marriott, ‘Tests of significance in canonical analysis,” 
Biometrika, v. 39, 1952, p. 58-64. 


Given two sets of variates x1, x» and (pb Sq) andu+1 
observations on each set, the canonical correlations are defined as the p roots 
(il? = 1? = --- 2 1,*) of the determinantal equation |Q — PT| = 0, where 
T is the dispersion matrix of the y’s; W is the dispersion matrix of the y’s 
with the x’s eliminated; Q = T — W is the dispersion matrix due to regres- 
sion. 

The exact null distribution of the largest root /;? is given by Roy! in a 
recursive form involving multiple integrals. The author gives the exact 
distribution of /,? for p = 2 and p = 3, g = 4. A significance test which is 
exact for practical purposes is given for p = 3and p = 4,g = 5.5% and1% 
significance levels to 2D (3 or 4S) for 3” /,? (m large) are given for p = 2, 
q = 2(1)12, 21; p = 3, g = 3(1)12, 21; p = = 5. 

An approximate test, 


Xion = — — + + 1)} log (1 — 1%) 


where D = p+q—1+4{(p — 1) (¢ — 1)}!, based on Wirks’ criterion? 
is proposed. 5% and 1% significance levels to 2D (3 or 4S) for $” 1,? (n large) 
are given for p = 2, g = 2, 6, 12, 21; p = 3,g = 3, 6,12, 21;p = 4,¢g=5 
which compare favorably with the exact values. The derivation of the test 
is omitted, with reference to the author’s unpublished thesis.* The 5% points 
of 1,? for the exact and approximate tests for p = 2, g = 5 are given to 2 or 
3D for m = 10, 20, 50, 100, ©, which indicate good results for 2 down to 20. 
I. OLKIN 

Michigan State College 
East Lansing, Michigan 

1S. N. Roy, “The individual sampling distribution of the maximum, the minimum, and 
any intermediate of the p-statistics on the null hypothesis,” Sankhyé, v. 7, 1945, p. 133-158. 


2S. S. Wixks, “Certain generalizations in the analysis of variance,” Biometrika, v. 24, 
1932, 471-494. 
3 


. H. C. Marriott, The Analysis and Interpretation of Multiple Measurements. 1951, 
University of Aberdeen. 


1054[K ].—Sice1t1 Moricut1, ‘‘A lower bound for a probability moment of 
any absolutely continuous distribution with finite variance,’’ Annals 
Math. Stat., v. 23, 1952, p. 286-289. 


The n-th probability moment of a population with probability density 
function f(x) is defined as 


a, « f dx. 


The author considers populations with zero mean and finite variance o?. A 

table is computed for the lower bound of the reduced probability moment 

Q,0° for m = 2(1)10 to 5D and a population which attains the lower 

bound is exhibited. The problem arose in the distribution of sample ranges. 
J. R. VATNSDAL 


State College of Washington 
Pullman, Washington 
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1055(K }.—Se1jt Naseya, “Absolute moments in 2-dimensional normal 
distribution,”’ Inst. Stat. Math., Annals, v. 3, 1951, p. 2-6. 


Let x and y have the joint probability density function 
1 x* 2pxy 


1 
220102 v1 = 2(i — \ 0102 


This paper presents expressions for E(|x"y"|) which are functions of o:, 
o2, and p. The cases considered are those for which simultaneously m = n, 
m2=1i1,2 = 0, m + 2 = 12. 


J. E. WaALsH 
U. S. Naval Ordnance Test Station 
China Lake, California 


1056(K ].—K. R. Narr, “Some three-replicate partially balanced designs,” 
Calcutta Stat. Assn., Bulletin, v. 4, 1951, p. 39-42. 


Four types of two-associate and four types of three-associate designs, 
using three replicates of each treatment, are discussed. A list of 48 useful 
designs is presented for which v (number of treatments) = 10 and & (number 
of treatments per block) = 10. The number of associate classes (classes of 
treatment comparisons) and the efficiency factor are given for each design. 

R. L. ANDERSON 
North Carolina State College 
Raleigh, N. C. 


1057[(K ].—M. H. QuENouILLE, “The variate-difference method in theory 
and practice,” International Stat. Inst., Review, v. 19, 1951, p. 121-129. 


The author discusses the usual use of the variate difference method to 
estimate the random variance in a time series and to test for the existence of 
trend after the i-th difference. He develops a new test of the existence of 
trend and a formula to estimate the random variance based on differences of 
the variances of the successive variate differences. If we let V; be the variance 
of the i-th set of variate differences (as defined by the author), then the 
variance of A"V; is amiot/n, where nm is the number of observations, o? is 
the true random variance, and values of a»; are given by the author for 
4 = 1(1)10 and m = 0(1)5, m + i = 10, to 4D for m = 0, 1 and to 6D for 
m = 2. 

The author often interchanges his notation (e.g., m is an order of differ- 
encing and also indicates both a terminal and an intermediate variance, 
V,,). He does not indicate how to decide on the best of many estimates to 
use for the random variance. However, his methods may be quite useful, 
if a systematic procedure is developed for using them. 

Some discussion is made of the difficulties involved when the residuals 
are serially correlated. The author discusses two empirical examples not 
involving serial correlation and one practical example with serially correlated 
residuals. More discussion is needed of how the empirical data were gener- 
ated. 

R. L. ANDERSON 
North Carolina State College 
Raleigh, N. C. 
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1058[K ].—P. V. Suxuatme, V. D. THawant, V. G. PENDHARKAR, & N. P. 
Natu, “Revised tables for the d-test of significance,’’ Indian Soc. of 
Agricultural Stat., Jn., v. 3, 1951, p. 9-23. 

The d-test under discussion is the Behrens-Fisher test for the difference 
of two sample means drawn from normal universes with possibly unequal 
variances. Specifically, 


d = (% — %,)/(s? + 51°)! = t2 cos — th sin 


in which Z, and Z2 are the sample means, s;? and s,* are unbiased estimates of 
the population variances, tan = s,/s2 and and are Student-Fisher 
with m, and m2 degrees of freedom. SUKHATME’s original tables of 5% and 1% 
values of d appeared in the 1943 edition of FisHer & YATEs’ tables.! Mean- 
while Fisher developed an asymptotic method of calculating these percent- 
age points? and detected a small positive bias in Sukhatme’s values. The 
latter investigated further and found that his errors were due to his use of 
linear interpolation in Student’s original tables of the é-distribution.* His 
corrected 5% values were used in the 1948 edition of the Fisher & Yates’ 
tables but his 1% values were corrected by Fisher for a positive bias shown 
in comparison with those of Fisher. In the present paper the authors report 
the results of a further investigation. They find that the positive bias re- 
ported by Fisher did not exceed unity in the third decimal place and that the 
use of quadratic interpolation is sufficient to give accuracy to 3D for the 
arguments used. Their investigation of Fisher’s method casts doubt on its 
adequacy at all points. Their finally corrected values appear as their Tables 
III and VII. These are for the same arguments as those in the Fisher & 
Yates tables and are also to 3D. (Values are apparently given to 4D when 
the final figure is 5 but no indication is given whether such a 5 is 5+ or 5-.) 
The 5% values agree with the 1948 Fisher & Yates values but there are 
three discrepancies in the 1% values as follows: 


Ne Sukhatme Fisher & Yates 
24 24 15°(75°) 2.784 2.785 
6 6 30°(60°) 3.558 3.557 
8 8 30°(60°) 3.241 3.239 


1R. A. FisHer & F. YATEs, Statistical Tables rag Biological, Agricultural and Medical 
Research. London. 2nd ed. 1943, 3rd ed. 1948. Table V3. 
2R. A. FisHer, “The asymptotic approach to Behrens’s integral, with further tables 
for the d-test of significance,” Ann. Eugenics, v. 11, 1941, p. 141-172. 
‘ . again “New tables for testing the significance of o! tions,” Metron, v. 5, 1925, 
p. 


1059[L ].—D. Ca ‘‘Complementi analitici e numerici allo studio delle 
aste vibranti,”” Accad. naz. dei Lincei, Rome, cl. sci. fis. mat. nat., 
Rendiconti, s. 7, v. 12, 1952, p. 76-83, 277-285. 

Let B,, p = 0, 1, --+ be the roots of cosh x cosx = 1, p = Oand p=1 
being the double root. Table I, p. 82, gives values of 8 ,,(—1)*[8,—(p—4) x], 
and sinh 8,/[cosh 8, + (— 1)*] for p = 2(1)5. 

Table V, p. 282, contains similar information for the roots of 


cosh x cosx = — 1. 
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Six other tables give values of certain integrals occurring in the physical 
problem. 
A. E. 


1060[L].—M. A. DeNnGLER, M. GoLanp & Y. L. LuKE, Tables of the func- 
tions f " (e« — 1) Vi + 6 du/u. Midwest Research Institute, Report 
August 1952, 19 p. 

This report tabulates 


A(y, 6) = — cos u) (u? + du 


By, b) = (a2 + 098 


for 
b = 0(.2)1.8, y = 0(.05)2(.1)2.5 


and 
b = 2(.2)4, y = 0(.05).5(.1)2. 


The introduction describes the computation. 9D were carried, and the 
results have been rounded to 7D for presentation in these tables. The entries 
should be correct to within one unit in the seventh place. Five-point Lagran- 
gean interpolation in the y direction gives results correct to about twounits 
in the sixth place. If b > 1, similar interpolation in the b-direction is reliable 
to about one unit in the fifth place. If 6 < 1, interpolation in the b direction 


is poor. 
A. E. 


1061[L].—C. Ferrari, “The turbulent boundary layer in a compressible 
fluid with positive pressure gradient,’ Jn. Aero. Sci., v. 18, 1951, p. 
460-477. 

There is a table (p. 476) of Ga(x) = 2?"n! i?" erfc x for x = 0(.05)1(.1)2.2, 

n = 2(1)5. The values are to 4D but “the fourth place is not accurate.” 

This table extends one by HARTREE.! 

D. H. L. 


1D. R. Hartree, “Tables of the error function,” Manchester Lit. and Phil. Soc. 
Memoirs and Proceedings, v. 80, 1935, p. 85. [Reprinted in CaRsLAwW & JAEGER, Conduction of 
Heat in Solids, p. 373, Appendix TIT} 


1062[L ].—H. J. Gopwtn, “‘A method for the evaluation of 
x™ exp (— a)" dx,” 


Quart. Jn. Mech. Appl. Math., v. 5, 1952, p. 109-115. 
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Put 


(2/x)' exp (— 44) dt = g(a), 


x™ {g(x)}"dx = Imn, 
P(m,0,x) = 
P(m,1r,x) = rxP(m,r — 1,x) + P’(m,r — 1, x), 
Ym,o = P(m, m + 20,0) + 20) ! 
The author computes ym,, and J,,, by expansion in inverse factorial 
series. He gives 9-10D tables of Im, for m = 0,1, 2, m = 1(1)20, and 9-15D 
tables of ym for m = 0, v = 0(1)18; m = 1, 2, v = 0(1)10. 
A. E. 


1063[L].—V. R. THIRUVENKATACHAR & B. S. RAMAKRISHNA, “A case of 
combined radial and axial heat flow in composite cylinders,’ Quart. 
Appl. Math., v. 10, 1952, p. 255-262. 


On p. 261 there is a 2D table for the first four roots of the equation 


Ji(x) J1(3y/2) Yily) — Jily) 
Tolx) Yo(y) — Joy) ¥u(3y/2) 


y = [10x? + = 1, 3,5. 


10x 


where 
A. E. 


1064[V ].—Rosert W. SmitH, Jr., HELEN E. Epwarps & Stuart R. 
BRINKLEY, JR., Tables of Velocity of Steady Laminar Flow in Channels of 
Rectangular Cross Section. Report of Investigations 4885, U. S. Dept. of 
the Interior, Bureau of Mines, July 1952, 41 p. 


Table 2 contains the solution of the boundary value problem: w22+ w73 
= — 2 in the rectangle R(|x| <1, |g| < 1/c), w = 0 on the boundary of 
R. For each o = .1(.1)1, w(x, y) is presented as a function of x = 0(.1).5(.05) 
.7(.02)1 and y = coy. In the first 20 columns y = 0(.1).5(.05).7(.02).9; in the 
remaining columns y = .905(.005)1 for o = .1 and .2, y = .91(.01)1 for 
o = .3(.1).7, and y = .92(.02)1 for o = .8(.1)1. 


1 1 
For the same range of o Table 1 has wo = w(0,0), 6 = f f w(x, y) 
0 0 


X dxdy, k’ = owo/8, k = o&/8; Table 3 has (— dw/dx).~1 for the above pairs 
of values of o and y; and Table 4 has (— dw/dy),-1 for the above values of x. 

In hydrodynamics w is the reduced velocity of the steady laminar flow 
of a viscous fluid through a channel with the rectangular cross section R. 

In the problem of twisting an isotropic prismatic bar with cross section R 
and shear modulus yu by a couple M whose moment is directed along the axis 
of the bar, w is the stress function. Moreover, if a is the angular twist per 
unit length of the bar, M = 8yuaa/o and the stress components are 7,2 
= cpadw/dy and = — padw/dx.! 

The tabulated functions were computed to 10D from Fourier series and 
rounded off to 6D. 
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The very last factor in equations (17) should be [g(0, y, «) P. The second 
of equations (12) is k’ = ---. The denominators of the factor before > in 
equations (5), (7), (9), and (14) are x°, and 

R. R. REYNOLDS 
NBSINA 


'1.S. SOKOLNIKOFF, Mathematical Theory of Elasticity. New York, McGraw-Hill, 1946. 

See under “Torsion” in index for references to electrical, hydrodynamic, and membrane 

ne ies. The right side of equation (38.15), p. 149, of the above reference is the series for 
x? +- 9*)/2 provided y, a, and 6 are replaced by 9, 2, and 2/c. 


MATHEMATICAL TABLES—ERRATA 


In this issue references have been made to RMT 1043, 1058, 1061, and 
1064. 


218.—N. W. McLacaian & P. HumBERT, Formulaire pour le Calcul Sym- 
bolique. Mémorial des Sciences Mathématiques, fasc. 100, 2nd ed., 
Paris 1950. 


Errata in the first edition have been noted in MTE 205 [MTAC, v. 6, 
p. 100-101]. Beside these errata the second edition contains also the follow- 
ing: 
. 4, formula 5; for #/4 read #/2 
. 4, formula 9; for E read F 
. 6, formula 8; for t* read t-* 
p read p 
a a 
. 12, formula 7; for — ¢’ (— log p) read ¢(log p) — f(0) 
. 12, formula 9; for 2" read 2? and for «/2Nt read x/V2t 
16, formula 9; add R(v) > 1 
21, formulas 6, 7; add 0 <t <1 
21, formula 11; add t >a 
. 22, formula 10; for ch?*t read sh*"*+'t (see Supplément, p. 10) 
. 22, formula 11; for sh?**"t read shvi 
. 22, last formula; delete the brackets [ ] 
. 23, formula 11; for p log (p/Vp* — 1) read p log (V p* — 1/p) 
25, formula 6; the I.h.s. is equal to » sh(v arg ch #) 
. 27, formula 6; for 1/(p + 1) read p/(p — 1) 
. 27, last formula but one; first term, for numerator on r. -h. s., read a 
. 28, formula 6; for R(v) > — 4 read R(v) > —1 
. 28, third last formula; for 2” read 2* and for R(u+¥v) >0 read 
Ru+v>-1 
. 28, last formula; for J,(t) read J2,(t), for (2vq — p) read (2vq + p) and 
for R(v) > 1 read R(v) > } 
. 29, formula 5; add R(v) > — 2 
. 29, formula 9; multiply the r.h.s. by 2! (see Supplément p. 11) 


. 11, formula 2; for 


. 23, formula 12; for p log 
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p. 30, third last formula; in numerator read (— 1)™a’t+?™hi@+)+m 
xr (Z+m+1) and add R(v) > — 2 

p. 31, formula 5; for J, read J,? 

p. 35, third last formula; for R(u + v) > 1 read R(u + v) > — 1 

p. 35, last formula; for I, read Ix), for (2vu — p) read (2vu + p), and for 
R(v) > 1 read R(v) > 4 

p. 36, formula 3; delete 

p. 36, formula 4; on r. » s. alter center sign to + and add R(v) > — 2 

p. 36, formula 10; for x I,42(2/p) read (—1)°J,4,(2/p) 

p. 36, last two formulae ; delete 

p. 37, formula 9; replace r.h.s. by pvp +a- Vp)”"/ va’, and add 
R(v) > 0 

p. 37, formula 10; for e—* read e~** 

p. 37, formula 12; for ap/s read ap/s* 

p. 39, formula 9; for } > R(v) > — } read — 1 < R(v) <1 

p. 39, formula 11; to l.h.s. add — 3 log {(¢ — b)/(¢ + 6)} Io(ay) 

p. 39, formula 12; to I.h.s. add — 5 log {(¢ — b)/(t + b)} Io(ay) and for 


pin[ ]read — p 
. 39, formula 13; delete 
. 39, formula 14; for I.h.s. read 


p. 40 fourth last formula: for — i read 2Vit and for (— )"/? read ei**/2 
p. 42, formula 3; to l.h.s. add — 3 log {(¢ — b)/(t + b)} (ber ay + 7 bei ay) 
p. 42, formula 4; ; for the |.h.s. read 

—b 
t+6 


x { er, ay + i kei, ay) 


— (bers ay + i bei-ray)| 1 


42, formula 5; to I.h.s. add — 5 (53) (ber ay + bei ay) 


. 42, formula 5; delete 

. 47, formula 4; for — a*/s* read — a?/s* 

48, last formula; for p?"*' read p*" 

50, third last formula; for read 

51, replace this page by p. 53 of the Supplément 

55, formula 3; for 2" read 2”! and for a/2Vit read a/V2t 
55, fifth last formula; replace r.h.s. by p[(1 — 1/p)™ — 1] 
. 56, formula 13; for Lon read Lin 

57, formulas 1, 2; delete 

57, last formula; for [1 + Bp]***—' read [1 + Bp]*te+ 
. 58, formula 6; for + arc tg p read — arc tg p 


~ 2 sin vr 
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p. 58, last formula but one; for e~*” read e~*? 

p- 59, formula 1; for (1 + h/p) read (1 + (hp)—) 

p. 59, third last formula; for 2A; read A; 

I am indebted to A. Erpféty1 for many of these corrections, some of 
which were communicated to him by O. VOELKER. 


N. W. 
Vizand & Co. 


51 Lincoln’s Inn Fields 
London W.C.2., England 


219.—NBSMTP., Tables of Fractional Powers. New York, 1946. 
Table 3, p. 34, for x © = 1.0678289226- - - 
read = 1.0678279226---. 


Mur an S. CorRINGTON 
RCA Victor 
Camden 2, N. J. 


220.—B. VAN DER Pot, “On the non-linear partial differential equation 
satisfied by the logarithm of the Jacobi theta-functions, with arithmetical 
applications, I,’’ Nederl. Akad. Wetensch., Proc., s.A., v. 54 [Indaga- 
tiones Math., v. 13}, 1951, p. 261-284. 


p. 281 for B2s = 336 87218 32202 92775 96104 01280 
read Bos = 436 56892 24858 87663 46104 01280 
B. VAN DER PoL 
12 Chemin Kreig 
Geneva, Switzerland 


UNPUBLISHED MATHEMATICAL TABLES 


151[F].—A. GLopEn, Factorisation of N*+ 1 for isolated values of N 
between 30000 and 40000, II. Two manuscript pages. Deposited in the 
UMT Fite. 


This constitutes an extension of UMT 144 [MTAC, v. 6, 1952, p. 102] 
and gives 50 new factorisations. 
A. GLODEN 
11 rue Jean Jaurés 
Luxembourg 


152[F].—A. GLopEN, Table of the Least Solution of the Congruence 2x* + 1 
= 0(mod p*) and Factorisation of the Corresponding Numbers 2x* + 1. 
Three manuscript pages. Deposited in the UMT FILE. . 


The prime ? is taken less than 1000. 
The largest number 2x? + 1 factored is 


2(380552)? + 1 = 3-11-883*?-11257. 

A. GLODEN 

11 rue Jean Jaurés 
Luxembourg 
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153[F].—A. GLopEN, Factorisation Table for the Numbers N* +- 1, N = 500. 
Six typewritten pages. Deposited in the UMT Fine. 
The table is an extension of CUNNINGHAM’s! table to NV < 200. Of its 500 
numbers 147 are completely factored. All unknown factors exceed 600000. 
A. GLODEN 


1A. J. C. Connincuam, Binomial Factorisations. V.6, London 1923, p. 140-141. 


154[F].—F. GRUENBERGER, Lists of Primes. Two sheets tabulated from 
punched cards. Deposited in the UMT FILE. 


The list of primes is extended from 50039981 to 50060033. There are 
1131 primes between these limits. This is a continuation of a list given in 
UMT 148 [MTAC, v. 6, p. 167]. 

F. GRUENBERGER 
Univ. of Wisconsin 
Madison, Wis. 


155[F].—R. J. Porter, Tables of Irregular Negative Determinants of ex- 
ponent 3n. Typewritten manuscript on deposit in the UMT FILe. 


The table gives the values of D < 50000 for which there is a determinant 
— D which is irregular with an exponent of irregularity which is divisible by 
3. [See Dicxson’s History! for definition of these terms. ] 

The table is arranged by thousands. There are 11, 17, 21, ---, 43 D’s in 
the first, second, ---, 50th thousand, a total of 1718 D’s altogether. Most of 
these have exponent 3. Only D = — 17561 has an exponent 6. Thirteen 
however have exponent 9. These are — D = 3299, 6075, 11907, 17739, 
23571, 24300, 27675, 29403, 33075, 35235, 41067, 46899, and 47628. All 
other D’s have exponent 3. 

The list was constructed by making extracts from some hundreds of the 
writer’s series of determinants of class-number 3k. To each determinant in 
these series belongs a class which has the property of duplicating into its 
own opposite; e.g., the determinant 21481 has a class (149, 71, 178) which 
duplicates into (26522, 8117, 2485) and thence by reduction to (2485, — 662, 
185), (185, — 78, 149) and (149, — 71, 178). 

These extracts are filed in numerical order with their corresponding A 
values (e.g., 149 in the above) and any determinants which have more than 
one entry of A values against them are irregular (exp. 3m). 

It is found, in practice, that to make extracts from the series for each 
block of 10,000 determinants takes approximately 40 hours’ work. 

R. J. PoRTER 
266 Pickering Road 
Hull, England 


1L. E. Dickson, History of the Theory of Numbers, v. 3, Washington 1927, New York, 
1934, Chap. 5. 
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Edited by the Staff of the Machine Development Laboratory of the National Bureau 
of Standards. Correspondence regarding the Section should be directed to Dr. E. W. 
Cannon, 415 South Building, National Bureau of Standards, Washington 25, D. C. 
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TECHNICAL DEVELOPMENTS 
THE USAF-FAIRCHILD SPECIALIZED DIGITAL COMPUTER 


The USAF-Fairchild Computer (or “SPEC,” for Special Purpose Elec- 
- tronic Computer) was designed and constructed at the NEPA Project' in 
Oak Ridge, Tennessee, to obtain solutions for a specific type of problems. 
These problems consisted of large systems of linear, algebraic equations 
involving up to 300 variables. 

The SPEC, shown in the frontispiece, is an electronic, fixed-decimal 
machine, with all numbers lying between — 1 and + 1, and has four-digit 
capacity, although it will be explained later that eight digits are carried in 
the accumulator. By specializing the purpose of the computer, construction 
time, development of special components and circuitry were minimized. It 
was also possible to simplify the problem insertion so that little skill is re- 
quired of the operating staff. 


CONTROL 
UNIT 
PRINTER 
LONG 
TAPE 4 
MULTIPLIER ACCUMULATOR 
| 
TAPE 
TEMP. 
STORAGE 
UNIT BLOCK DIAGRAM OF COMPUTER 


Fic. 1. 


The NEPA Computer employs the Gauss-Seidel iteration method for the 
solution of simultaneous, linear equations. By this method, a sequence of 
approximate solutions is obtained by successive iterations such that the 
limits of these solutions are the correct solutions of the equations. 

When the problems have been arranged for solution, each of the equa- 
tions in the set is solved for one of the variables as a function of a constant 
term and the remaining variables. In this manner, if the order of the set is NV, 
there will be N equations of the following type: 


j=N 
(1) Xi =F AgX; + where Ay = 0. 


The method requires the summing of products and the addition of a 
constant term C;. Values are assumed for each variable and substituted 
successively in the equations where at the conclusion of the evaluation of the 
right-hand member of each equation the new value for X; replaces the old. 

In the SPEC, two synchronous, magnetic tapes provide the required 
storage. These are shown in the block diagram of Fig. 1. A short, endless 
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tape stores the unknowns. As this tape revolves, the unknowns are read in 
order. Values for the unknowns, altered by successive iterations, will be 
repeated during succeeding tape revolutions. All the coefficients and con- 
stants are stored on the long tape, with accompanying instructions to the 
control unit. An instruction indicates to the computer which multiplier is to 
be selected for the number being read from the long tape and also signals 
the conclusion of each equation, at which time the new result is printed. 

A common sprocket drive provides simultaneous reading of the two 
tapes. The synchronized drive permits arrangement of the coefficients on the 
long tape and the unknowns on the short tape in such a manner that a co- 
efficient and the unknown by which it is to be multiplied are read by the two 
inputs simultaneously when they are needed during the course of the com- 
putation. 

The computer is, basically, a summing device for a series of products. 
As shown in Fig. 1, the arithmetic unit consists of a multiplier coupled to an 
accumulator. The numbers received at the two inputs to the multiplier are 
multiplied and then added in the accumulator to the results obtained from 
previous similar operations. The control unit determines the sign of the 
operation to be performed by the accumulator, the product being either 
added to or subtracted from the total in the accumulator as the computation 
requires. 

The multiplier is a diode matrix incorporating in its connections all the 
products of two decimal digits between 0 and 9. A single digit from each 
input number is coupled to the matrix, and the output of the matrix repre- 
sents the product of those digits selected at the input. A complete product 
of two four-digit numbers requires sixteen operations, with the partiai prod- 
ucts being added, as obtained, into appropriate locations in the accumu- 
lator. The complete product requires approximately 400 microseconds. 

The accumulator capacity is eight decimal digits. The eight-decimal- 
digit capacity prevents the formation of an excessive error due to repeated 
roundoffs in the summation of a series of products. Roundoff of the result is 
accomplished only at the conclusion of a series, prior to the transfer of a 
result out of the accumulator. 

Each decimal digit is represented by a coded group of four binary digits. 
All operations in the accumulator are binary operations, using the ‘‘excess- 
three” code. The accumulator consists of 32 parallel-fed binary adder 
circuits. Sufficient inputs are available to handle the number and the 
required ‘‘excess-three’’ correction factor simultaneously. 

At the conclusion of each series of products, and upon command from 
the control unit, the resulting total is rounded off and transferred both to a 
printer circuit, which provides a typed copy of the results for the operator, 
and to the temporary-storage register. This transfer occurs at the end of a 
sequence of operations, and the short tape is not then in the proper position 
for recording the new X over the old X. The temporary-storage register 
holds this X value until the point on the short tape is reached at which this 
value is to be recorded. The recording will normally occur during the solution 
of the next equation, and the computer will both record and use the X in the 
temporary-storage register during the solution of that equation. The re- 
cording location on the short tape is determined by the set of instructions 
placed on the long tape by the operator prior to the solution of the problem. 
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A total of 2206 vacuum tubes are used in the computer. Over half of 
these tubes are high-vacuum diodes, types 6AL5, which were used in pref- 
erence to crystal diodes for matrix circuits. 

Design and construction of the SPEC started in July 1949, and the first 
problem was solved on February 1, 1950. It was in continuous operation at 
the NEPA Project from the spring of 1950 until February 1951 and is being 
installed at the Oak Ridge National Laboratory. The computer will be called 
the “ORACLE” which is derived from Oak Ridge Automatic Computer for 
Linear Equations. 

During the operating period, the computer was in an operating condition 
for approximately 85% of the total available time, with 15% of the time 
devoted to testing and servicing. Most machine errors were detected during 
the solution of problems, although test problems were used for periodic 
checks to assure proper operation of the machine. Check circuits, with visual 
and audible indications, are employed to detect improper operation of the 
machine in the most critical places, but no compiete checking system is 
employed. 

Throughout the period of operation, the computer was exceedingly use- 
ful as a means of obtaining solutions to systems of simultaneous, linear 
equations up to the limit of its useful capacity. In addition, problems in- 
volving matrix products, Fourier analysis, numerical integration, and 
matrix inversion were undertaken with considerable savings in time and 
effort. 

The author wishes to acknowledge the outstanding contribution of Mr. 
V. G. Lewis and his staff of technicians at the NEPA facility in Oak Ridge 
in the construction and assembly of the computer. Special recognition is due 
Mr. L. C. Oakes for his valuable contributions in testing, servicing and 
operating the machine. 

J. J. Stone 
Oak Ridge National Laboratory 
Oak Ridge, Tenn. 


1 Nuclear Energy for the Propulsion of Aircraft Project, conducted by Fairchild Engine 
and Airplane Corporation under special contract with the United States Air Force. 


DISCUSSIONS 


FLOATING OPERATIONS ON THE EDSAC 


Summary. The difficulties which arise when programming calculations 
for large automatic calculating machines which have a fixed decimal point 
are discussed. This leads to a consideration of the possibility of using 
floating decimal arithmetic for certain kinds of calculations. A method by 
which floating decimal arithmetic can be carried out with any fixed decimal- 
point machine is outlined and the scheme adopted for use with the EDSAC 
is described in detail. 

This scheme is based on a special kind of subroutine which we call 
interpretive. This enables the programmer to use a new order code of his own 
choice. The ‘orders’ of a programme drawn up with such a code are selected 
under the action of the interpretive subroutine and interpreted in terms of 
sequences of orders which perform the required operations. With the EDSAC 
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a single address ‘order’ code has been adopted, similar to the ordinary order 
code, in which the arithmetical ‘orders’ are interpreted in terms of floating 
decimal arithmetic. Special ‘orders’ for simplifying counting operations and 
the modification of other ‘orders’ are provided. The ‘order’ code is described 
in detail and an example is given of a programme drawn up using this code. 
Other topics discussed are the use of auxiliary subroutines with the inter- 
pretive subroutine, a method of facilitating programme assembly, and tech- 
niques for using the input tape as a form of auxiliary store. Finally the times 
of operation of the individual ‘orders’ of the code are given, together with an 
estimate of the factor by which the calculation time is increased as a result 
of using floating decimal arithmetic for an entire calculation. 

Introduction..One of the most tedious tasks arising in the programming of 
calculations for solution by fixed decimal-point machinery is to arrange the 
calculation so that all quantities concerned remain within the limits of the 
machine and yet are expressed to an accuracy sufficient to ensure the desired 
precision in the results. Our experience with the Epsac has shown that if 
this task can, in some way, be relegated to the machine then the time taken 
to prepare a programme after a problem has been understood is considerably 
redu 

The problem does not arise with machines designed to operate directly 
with numbers expressed in the floating radix form. Numbers in this form are 
represented by a-r?. The first machine of this kind was the Bell Telephone 
Laboratories Relay Computer Model V™. This is a decimal machine (that 
is,r = 10) in which 1 > |a| > 0.1,19 > » > — 19 and ais expressed to an 
accuracy of seven significant figures. Since this was completed all important 
relay machines have been equipped with similar facilities. No electronic 
machine of this kind has yet been built but we would remark that in our 
opinion an electronic machine provided with a floating point arithmetical 
unit would be a powerful computing instrument even if it had a relatively 
slow store, a magnetic drum, for example. It would be particularly suitable 
for a laboratory which had to solve quickly a wide variety of problems as 
they are presented. 

For a fixed decimal-point machine the usual arrangement is to associate 
scale factors with some or all of the quantities occurring in the calculation. 
These scale factors are of two kinds. First certain quantities can be asso- 
ciated with fixed scale factors which remain unaltered throughout the cal- 
culation. Their value must be chosen so that the quantities concerned do 
not exceed capacity and yet can be represented at all times with sufficient 
accuracy. This is the problem of providing ‘elbow room.’ It may not be pos- 
sible to satisfy both requirements by the use of a single fixed scale factor and 
it is then necessary to introduce an adjustable scale factor, that is, a scale 
factor which is altered during the course of the calculation. It may be 
adjusted either continually or occasionally in accordance with some criterion 
or it may take a preassigned set of values. All these things the coder has to 
think about if the calculation is to be arranged to use the least possible ma- 
chine time. The difficulties arising will vary from one calculation to another 
and may be trivial, moderately complex, or really hard. The individual na- 
ture of calculations and the necessity for a mathematical understanding of 
them, however, makes a general solution to the problem on these lines un- 
attractive to attempt. The only way in which the difficulties may be avoided 
is to associate every number occurring in the calculation with its own ad- 
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justable scale factor. The scale factors can be stored most economically if 
they are powers of 2 or of 10. In the latter case the above scheme virtually 
amounts to programming a floating decimal point. This is the solution that 
we feel should be adopted for many kinds of calculation. 

The floating decimal form of representation of numbers is especially 
useful when it is required to evaluate algebraic expressions. For example, to 
code the evaluation (au? + bu + c)/(du + e) with a fixed decimal-point 
machine can be very troublesome if it is required to maintain accuracy over 
the entire range of numbers. If the reader has any doubts about this he is 
advised to try it. If floating decimal arithmetic is used, then the coding 
requires little thought and is quite direct. 

It is the purpose of this paper to describe a method by which floating 
decimal arithmetic can be carried out on fixed decimal-point machines, and 
in particular to describe the scheme adopted for the EDSAC. 

General considerations. The following remarks are fairly general and apply 
to almost any machine although certain features mentioned, for example, 
short and long locations, are made with the EDSAC in mind. 

The most convenient way to programme operations on numbers ex- 
pressed in the floating decimal form is by means of a special type of sub- 
routine which we call interpretive. This enables words similar to those repre- 
senting ordinary machine orders to be interpreted in terms of floating decimal 
operations. Such words will be referred to as ‘orders’ (with quotation marks). 
Each ‘order’ resembles an ordinary machine order but is never obeyed as 
such by the control circuits of the machine. Instead, ‘orders’ are selected in a 
definite sequence by the action of the interpretive subroutine and interpreted 
in accordance with a preassigned ‘order’ code, by means of sequences of orders 
which form part of the interpretive subroutine. The selective action of the 
interpretive subroutine will be referred to as the ‘control’ (in quotation 
marks). The advantage of such interpretive subroutines is that the 
‘order’ code may be chosen to suit the convenience of the programmer. 
There need be no relation between the form of the ‘order’ code and 
the ordinary code of the machine. However, we have become so familiar 
with the ordinary EDSAC order code that a similar single address ‘order’ 
code has been adopted for the interpretive subroutines. Interpretive sub- 
routines have also been placed in the EDSAC library for carrying out 
arithmetical operations on complex numbers and on double-precision 
numbers. 

When representing floating decimal numbers in a fixed decimal-point 
machine it is most economical to pack both the numerical part and the 
exponent into a single storage location. This means, of course, that the 
numerical part of the number has fewer binary digits to represent it than 
would normally be available. However, they are all significant whereas in 
fixed decimal-point working digits are thrown away to provide ‘elbow room.’ 
In both cases accuracy may be lost owing to cancellation of leading digits at 
one end of a number and to rounding-off errors at the other. 

Two long and two short storage locations are set aside to form a kind of 
‘arithmetical unit.’ One long location holds the numerical part of a number 
and one short location holds the exponent. Together they form the floating 
decimal accumulator. In a similar fashion the other long location and the 
other short location form the floating decimal register. 
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An arithmetical ‘order’ first causes a subroutine to select and unpack the 
operand and place the numerical part and the exponent into a further pair of 
storage locations—the ‘multiplier register.’ In the case of an add ‘order’ the 
operand is then added to the number held in the floating decimal accumula- 
tor. This is done as follows. Let ao, a2, and a, denote the numerical parts of the 
operand, augend, and sum respectively, and let po, a, and p, similarly 
denote their exponents. Then we have 


a,107* = [ap + 070 => Pa 
= [ao10- + po < Pa 


The subtract ‘order’ works in a similar fashion. A register ‘order’ places the 
operand in the floating decimal register. Other ‘orders’ enable the product of 
the operand with a number held in the floating decimal register to be added 
to or subtracted from the number held in the floating decimal accumulator. 
The transfer order causes the number held in the floating decimal accumula- 
tor to be converted to standard form, packed, and finally transferred to the 
store; the floating decimal accumulator is then ‘cleared’ by replacing the 
number held in it by zero, that is, by the special number 010-*. The input 
‘order’ causes the two parts of a number to be read from the input tape, 
packed, and transferred to the store. The output ‘order’ causes the numerical 
part and the exponent of the number held in the floating decimal accumula- 
tor to be printed on the same line in two adjacent columns on the page of the 
teleprinter. 

The use of two separate storage locations for the floating decimal accumu- 
lator allows the range and accuracy of numbers held therein to be greater 
than those held in a single storage location elsewhere. This enables products 
to be accumulated without loss of accuracy due to intermediate rounding- 
off errors. 

The scheme adopted for the EDSAC. The interpretive subroutine that has 
been developed to carry out floating decimal arithmetic with the Epsac will 
now be described in more detail. Throughout this and the following section 
references to specific features or conventions used with the Epsac have been 
avoided as far as possible but the following terms are used at least once: 
initial orders, control combinations, preset parameters (not to be confused 
with the current parameter defined below), short and long locations and the E 
order. Detailed descriptions of all these features can be found in reference 2. 

The following abbreviations will be adopted : 


F(A) denotes the floating decimal number held in the floating decimal 
accumulator 

F(R) denotes the floating decimal number held in the floating decimal 
register 

S(mD)_ denotes the long storage location having address m 

S(mF) denotes the short storage location having address m 

F(mD) _ denotes the floating decimal number held in S(mD) 

F(mF) denotes the floating decimal number held in S(mF) 

C(mD) denotes the ordinary number held in S(mD) 

C(mF) denotes the ordinary number held in S(mF) 


Number representation. The method described above of packing the nu- 
merical part a and the exponent ? of a floating decimal number a-10? in a 
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single storage location has been adopted. With the Epsac it amounts to 
representing the floating decimal number by the ordinary (fractional) 
number p2-* + a2-’, where a is expressed with an accuracy of 28 binary 
digits and lies in the range (1 > |a| > 0.1) and # is an integer such that 
63 > p > — 63. Zero has the special representation 010-*. When numbers 
are transferred from the floating decimal accumulator to the store they are 
automatically transferred to the standard representation. 

The ‘order’ code. Since the ‘orders’ are similar in form to the ordinary 
machine orders of the EDSAC they will be presented in the conventional form 
adopted for the latter. The ‘orders’ of the code fall into two classes, arith- 
metical and organizational. The arithmetical ‘orders’ can refer either to short 
or to long storage locations according as the ‘order’ is terminated with the 
code letter F or D. This is the usual EDSAC convention. If, however, these 
‘orders’ are terminated by A instead of D or A instead of F then a number, 
known as the current parameter, held in a certain storage location, will be 
added to the address of the order before it is obeyed. Among the organiza- 
tional ‘orders’ are two—the P and the F orders—which enable the current 
parameter to be set to some initial value and subsequently adjusted after 
each cycle of a repetitive operation. In this way an arithmetical ‘order’ can 
refer to different numbers during different cycles of the calculation. These 
orders are described in more detail below. The scheme is similar in some 
respects to the way in which orders are modified in the Manchester Univer- 
sity Electronic Computer Mk. 2. In this machine a built-in facility enables 
a number held in one of 8 registers—called the B registers—to be added to 
the address of an order immediately before that order is obeyed. 


The arithmetical ‘orders’ 


AmD_ Add F(mD) to F(A) 

BmD_ Subtract F(mD) from F(A) 

HmD Replace F(R) by F(mD) 

VmD_ Add the product F(R)-F(mD) to F(A) 

NmD _ Subtract the product F(R)-F(mD) from F(A) 

DmD_ Replace F(A) by F(A)/F(mD) 

@mD_ Replace F(mD) by F(A) and F(A) by 0.10-* 

LmD_ Similar to ¢ m D but in addition the non-standard content of the 
floating decimal accumulator is printed as negative sign or space ; 
exponent (2 figs.); space; negative sign or space; 8 decimal 
digit fraction. For example, — 98.283742 would be printed as 
02 — 98283742. 

AmD Input a sequence of numbers on the tape terminated by X into the 
locations mD, (m — 2)D, etc. Each number is punched in the 
following way: Characters to represent the exponent; sign; 
numerical part. In the numerical part the decimal point is under- 
stood to be before the first digit punched. Any number of digits 
may be punched. For example, — 98.283742 would be punched 
as 2 — 98283742. X is punched before the first number of the 
sequence and then the sequence is copied on to the data tape in 
the reverse direction. Numbers read into the machine must have 
exponents in the range — 16 to + 15, so that it can be repre- 
sented by a single tape character. 
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The organizational ‘orders’ 


MmF _ Conditional order: if F(A) > 0 transfer ‘control’ to the ‘order’ 
which stands in S(mF) ; otherwise replace F(A) by | F(A)| and 
proceed with the next ‘order.’ 

Cm F Transfer control of the machine to the order which stands in 
S(mF), that is, return to the machine order code. 

XmF _ _ Transfer ‘control’ to the X-auxiliary whose entry ‘order’ stands 
in S(mF). When the X-auxiliary is completed it will return 
‘control’ to the ‘order’ immediately following the X ‘order.’ 

GmF Transfer ‘control’ to the ‘order’ standing in S(mF). 

a ; The function of these ‘orders’ is described below 

The P and F ‘orders.’ These orders are provided to facilitate the coding of 
cycles of orders and of cycles within cycles, etc. Each cycle of ‘orders’ or 

‘loop’ is started by a P ‘order’ and terminated by an F ‘order.’ The two are 

said to correspond. This correspondence is similar to that between left and 

right handed brackets in a lengthy algebraic formula if the direction left to 
right in the formula corresponds to the direction in which ‘control’ is nor- 
mally advanced. Associated with each loop of ‘orders’ is a parameter whose 

value is set, initially, by the P ‘order’ and is subsequently adjusted by the F 

‘order.’ The current parameter (see above) associated with any particular 

‘order’ of a loop is the parameter set by the P ‘order’ which the ‘control’ 

encountered last. It is now possible to give a formal description of the func- 

tion of these ‘orders.’ 


PmF Record a new parameter—the current parameter of the following 
‘orders’—and set it to the value — m2-% 

Fm F If thecurrent parameter is < — m2-, increase its value by m2-%, 
return ‘control’ to the ‘order’ immediately following the corre- 
sponding P ‘order’; otherwise proceed with the following ‘order’ 
after restoring the current parameters to the value it had before 
the corresponding P ‘order.’ 


The above description should be sufficient to enable a programmer to 
use the P and F orders in a programme. We now give a detailed description 
of the means by which the effects of these ‘orders’ are achieved. 

The number of parameters has been limited to 3. This restriction, al- 
though not essential, is reasonable as problems involving loops within loops 
more than three ‘deep’ are likely to occur very rarely. 

The parameters are stored together with the addresses of the locations of 
the corresponding P ‘orders,’ in 3 long storage locations S(hD), S(h + 2)D, 
and S(h + 4)D. At any stage in the programme S(hD) contains the current 
parameter, S(h + 2)D the previous parameter, and S(h + 4)D the previous 
parameter but one. 

A P ‘order,’ for example, P m Fin S(nF) causes C(h + 4)D to be replaced 
by C(h + 2)D, C(h + 2)D to be replaced by C(hD), and C(hD) to be re- 
placed by the number — m2- + n2-**. When the corresponding F ‘order,’ 
F q F is encountered the following operations take place. The size of the 
current parameter (the number — m2-* in the short storage location 
S(h + 1)F) is tested. If this number is less than — g-2-" it is increased by 
the amount g-2-"* and ‘control’ is transferred to the ‘order’ which stands in 
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S(n + 1)F, where n-2-" is the number in the short location S(hF). Other- 
wise C(4D) is replaced by C(h + 2)D, C(h + 2)D is replaced by C(h + 4)D 
and then ‘control’ proceeds with the next ‘order.’ 

Example. The action of the ‘orders’ and the context in which they are 
meant to be used may be more clearly understood from a study of the follow- 
ing programme for the calculation of a root of an algebraic equation by the 
Newton-Raphson iterative process. 


Location of data: \et f(x) = Er a,x"-* = 0 be the equation to be solved. 
r=9 


It is assumed that the coefficients dn, @n-1, ---, @ are in S(100D), S(98D), 
-++, 5(100 — 2n)D. The initial approximation x» stands in S(6D) and all 
subsequent approximations are placed there. A small quantity 5, used in the 
convergence criterion, stands in S(8D). 

Formula used : the iterative formula 


= Xn — f(%n)/f (xn) 
is used. The iteration is arranged to terminate when |x, — x»1| < 4. 
f(x) is calculated from the recurrence relations 
Go = Go; = + In = f(x) 
and f’(x) is calculated from the similar recurrence relations 


go = 0; = Gr + Gr; gn = f'(x). 
The orders of the programme are listed below. For those readers who are 
not familiar with the Epsac conventions it may be stated that the code 
letter 0 provides for a system of numbering relative to the order immediately 
following a control combination G K. Thus, in the example below, the ‘order’ 
M 18 @ refers to the ‘order’ B 8 D. 


G K 
0 ¢ D ‘clear’ floating decimal accumulator 
1 © 10 D ‘clear’ ; S(10D) 
2-H 6 D place x, in the floating decimal register 
3 P 2n4+2 F 
4 12 D 
5 V 10 D 
brit’ = + cycle n + 1 times 
10 F 2 F 
10 Dt Ges) to 
13 B 10 D 
14 A 6 D = Xn — f(xn)/f' (xn) 
15 @ 6 D 
16 A 10 D oO 
17 M 18 6 form modulus of (xn41 — a) 
18 B 8 D subtract 6 
19 M 1 6 discriminate 
20 @ D ‘clear’ floating decimal accumulator 
~4 4 } print the root x 
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Auxiliary subroutines. To facilitate the coding of entire calculations in 
floating decimal arithmetic it is useful to extend the range of action of 
the interpretive subroutines by the addition of auxiliary subroutines for 
carrying out common numerical processes. The auxiliary subroutines, hence- 
forth referred to simply as auxiliaries, are of two kinds. The first kind con- 
sists entirely of ordinary orders and will be referred to as C auxiliaries. 
Auxiliaries of the second kind consist largely of ‘orders’ and hence themselves 
use the interpretive subroutine. They will be referred to as X auxiliaries. 

C auxiliaries. These are called in by a C ‘order’ which transfers control to 
the entry order of the auxiliary. This then carries out the appropriate calcula- 
tion by means of machine orders. When this has been completed control is 
transferred back to a point in the interpretive subroutine. This causes the 
‘control’ to resume the obeying of ‘orders’ starting at the ‘order’ immediately 
following the C ‘order.’ Thus the C auxiliary can be considered as being an 
extension of the ‘order’ code and used as such. 

A typical C auxiliary is one which replaces F(A) by its square root. The 
way in which this is done using ordinary orders is as follows. Let p, and mm. 
denote the exponent and numerical part of F(A). Two cases arise: If p, is 
even, p, = p./2 and n, = Va; if pa is odd, then p, = (p. + 1)/2 and n, 
= ¥n,/10. In both cases the arithmetical operations can be done most simply 
by using the ordinary techniques of fixed decimal-point working. 

C auxiliaries are almost as fast and as economical in storage space as the 
corresponding routines for fixed decimal-point working. The extra orders 
required to handle the exponent are largely offset by the orders that would 
otherwise be required to cater for a large range of the numerical part. 

X auxiliaries. These are called in by the special ‘order’ X m F, where m 
is the location of the entry order of the auxiliary. When the X auxiliary has 
completed its part of the calculation ‘control’ is returned to the ‘order’ 
immediately following the X ‘order’ which called the auxiliary into use. This 
cannot be done by the same method as is used for the C auxiliary because the 
X auxiliary itself uses the interpretive subroutine. Instead the X ‘order’ 
causes the address of the location of the next ‘order’ to be stored in a certain 
location for future reference. This record is called the link. At the end of the 
auxiliary is a C ‘order’ which directs control to a set of orders within the 
interpretive subroutine which use the link to return ‘control’ to the main 
programme. 

X auxiliaries can use other X auxiliaries. To enable this to be done the 
X ‘order’ causes a list of links to be kept. At any point in the programme the 
link at the head of the list refers to the ‘order’ to which ‘control’ must be 
returned after the current subroutine has been finished with. The number of 
links in the list measures the depth to which ‘control’ has passed within the 
auxiliaries. A reference number which is adjusted every time a link is added 
to or removed from the list records this depth. This enables the link at the 
head of the list to be selected by those orders of the interpretive subroutine 
which are called into use when the X auxiliary has been completed. 

The Directory. Reference to the auxiliaries can be facilitated, if desired, 
by using the following scheme. The auxiliaries are first enumerated in the 
order in which they are to be read from the input tape into the store. Thus 
auxiliary no. 1 is the first auxiliary read from the tape, auxiliary no. 2 is the 
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second auxiliary read from the tape, and so on. When the auxiliaries have 
been ordered in this way the m-th auxiliary can be called in by the ‘order’ 
X m L or Cm L (depending on whether it is an X or a C auxiliary). 

This is achieved by means of a table of switching orders—called a 
directory—stored in consecutive storage locations beginning with S(hF). 
Each entry directs ‘control’ (or control) to the first ‘order’ (or order) of one 
of the auxiliaries. If, for example, the m-th auxiliary is an X auxiliary whose 
first ‘order’ stands in S(nF), then the m-th entry in the directory, that is, 
the entry standing in S(h + m)F, is the ‘order’ gm F. Thus to call in the 
auxiliary ‘control’ is first transferred, by the ‘order’ X h + m F, to the ‘order’ 
Gn F which in turn transfers ‘control’ to the first ‘order’ of the auxiliary. 
Similarly, if the m-th auxiliary is a C auxiliary, the m-th entry is the ordinary 
order E n F, and the auxiliary is called in by the ‘order’ Ch + m F. 

The directory is assembled by the initial orders as the auxiliaries are read 
from the tape into the store. At the same time the address of the locations, 
S(hF) of the first entry in the table is placed in the preset parameter location 
L so that X m + h F and Cm + h F can be punched as X m L and CmL 
respectively. 

The advantage of the scheme is that the master routine can be drawn up 
in the final form once the set of auxiliaries have been ordered. It is not 
necessary for the coder to keep a record of the locations in the store of indivi- 
dual auxiliaries. In effect, the scheme shifts the burden of ‘book-keeping’ 
from the programmer to the machine. Similar principles are used by the 
assembly subroutines (for programme assembly) which already exist in the 
Epsac library of subroutines (see reference 2). 

The use of the input tape as a form of auxiliary store. The Epsac is not 
provided with an auxiliary store so that for many problems the shortage of 
storage space is a real difficulty. This difficulty may be partly overcome by 
using the paper tape input medium as a form of auxiliary store. There are 
two ways in which this may be done. 

The first way, used when the storage requirements grossly exceed those 
available, is to carry out the calculation in stages, using the input tape to 
store numbers and ‘orders’ not required throughout the calculation. This 
method will be referred to as piecewise control of the calculation. 

In the second method the ‘orders’ of the master routine are placed on the 
tape and read into the store one at a time, each being obeyed immediately 
after it has been read. This method is used when the master routine consists 
of a lengthy sequence of ‘orders’ not many of which are repeated. If none of 
these are repeated the scheme takes no longer than putting the complete 
master routine into the store and then entering it. The principal advantage 
is that no storage space need be allocated to the master routine. A further 
advantage is that the progress of the calculation is apparent from the prog- 
ress of the tape. This method of control will be referred to as input control. 

In the Manchester University Computer Group, tapes drawn up on these 
lines are referred to as job-steering tapes. 

Piecewise control. To facilitate this mode of working the A ‘order’ and the 
initial orders can be used as follows. 

The A ‘order’ enables sequences of numbers of any length to be read 
from the tape into the store when required, overwriting, if necessary, informa- 
tion no longer wanted. 
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If the initial orders are retained intact during the course of the calcula- 
tion they can be recalled into use by a C ‘order’. In this way ‘orders’ (or 
orders) can be read from the input tape into the store in the usual way at 
any time during the calculation. The reading of ‘orders’ can be halted and 
control transferred back to the interpretive subroutine by a suitable control 
combination. 

The above scheme together with the use of the A ‘order’ enables quite 
complicated problems to be tackled. 

Input conirel. In this scheme a special C auxiliary—the input control 
auxiliary—is used. This causes ‘orders’ to be read from the tape and inter- 
preted immediately after they are read. Normally they will be of an arith- 
metical character or will call in auxiliary subroutines. The C auxiliary will 
continue to read and obey such ‘orders’ until this mode of ‘control’ is ter- 
minated by a suitable C ‘order’ on the tape itself. This directs control to a 
point within the interpretive subroutine and in this way causes ‘control’ to 
be returned to the ‘order’ immediately following the C ‘order’ which called 
the input control auxiliary into use. 

Times of operations of the ‘orders.’ The times of execution of the individual 
‘orders’ are as follows: 


‘orders’ Times of execution 
A,B 90 ms. 
V,N 105 ms. 
6G, MP. F 50 ms. 
D 140 ms. 
? 80 ms. 


The times of operation of the Z and A ‘orders’ are largely determined by 
the speed of the input and output units. The teleprinter can print six 
characters per second and the tape-reader can read about 25 characters per 
second. These rates allow for the time taken for binary to decimal conversion 
and vice versa. 

For comparison it should be stated that each ordinary order of the Epsac 
takes 13 ms. with the exception of the multiplication orders, V and N, which 
take 6 ms. each. 

From a direct comparison it would seem that the floating ‘orders,’ other 
than those used for reading and writing, are about 60 times as slow as the 
machine orders and hence that a programme using the interpretive sub- 
routine would be slower by the same factor. This is not altogether true be- 
cause in such a programme fewer ‘orders’ are needed than would otherwise be 
necessary as there are no scale factors to deal with and the techniques for 
counting and for the modification of ‘orders’ have been streamlined. More- 
over, the time taken by the C auxiliaries is about the same as that taken by 
the corresponding subroutine in fixed decimal-point working. 

These factors vary from problem to problem but our experience has 
shown that the reduction in speed varies from about 20 to 1 to about 4 to 1. 
The reduction of the time taken to code a problem has to be experienced to be 
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BIBLIOGRAPHY OF CODING PROCEDURE 


The National Bureau of Standards is forming a central library of notes, 
reports and technical publications concerned with programming and coding 
for electronic digital computers. 

The collection of material has in view the possibility of increasing the 
efficiency in the use of high speed digital computers. The material will include 
for example coding manuals, computing routines, supervisory routines and 
codes for generating other codes. Those interested in research and instruction 
in coding practices are invited to avail themselves of this material. 

Computation laboratories are invited to submit material to this library. 
A list of such acquisitions together with a short descriptive review will appear 
in the future issues of MT-AC; they will be numbered serially for reference. 
Material should be sent to the attention of J. H. WeGstEIn, Computation 
Laboratory, National Bureau of Standards, Washington 25, D. C. 

Some material of this kind has already been included or reviewed in 
MTAC. This bibliography begins with references to eleven such items. The 
last four items are new. 
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ting machine,” Jn. Sci. Inst. and Phys. in Industry, v. 26, 1949, p. 217- 
220. [MTAC, v. 4, p. 116] 

5. ANON., Description and Use of the ENIAC Converter Code. Ballistic 
Research Laboratories, Technical Note no. 141, Aberdeen Proving 
Ground, Maryland, November 1949, 23 pages, mimeographed. [MTAC, 
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6. R. W. Hammine, “Error detecting and error correcting codes,” Bell 
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9. M. V. WiLkEs, D. J. WHEELER & S. GILL, The Preparation of Programs 
for an Electronic Digital Computer. Addison-Wesley Press, Inc., Cam- 
bridge, Mass., 1951, 167 p., 22.9 K 15.2 cm. Price $6.00 [MTAC, v. 6, 
p. 51] 

10. R. A. Brooker, S. Git & D. J. WHEELER, ‘“‘The adventures of a 
blunder,” [MTAC, v. 6, p. 112-113] 

11. Davin P. Perry, “Minimum access programming,” [MTAC, v. 6, p. 
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UNIVERSITY OF ILLINOIS RESEARCH BOARD, University of Illinois 
Electronic Digital Computer Library Routines. This volume is an ex- 
pandable library of the routines which have been prepared for the 
Illinois machine and at present contains perhaps fifty routines. It 
contains among many other such routines as the Wheeler Leapfrog Test, 
the solution of simultaneous linear equations, solution of a system of 
first order differential equations, post mortem, integration routine, ad- 
dress search routine, eigenvalues and eigenvectors of a symmetric matrix, 
and square root. 

13. ECKERT-MAUCHLY COMPUTER CORPORATION, Programming Manual for 
the UNIVAC System. This describes the UNIVAC and how to code for 
it. The manual includes chapters on : flow charts, an illustrative problem 
(Denominated Payroll routines), floating decimal and double precision 
routines, elementary mathematical problems, input and output opera- 
tions, editing, problem preparation. 
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1000. MarTIN GARDNER, “Logic machines”, Scientific American, v. 186, 
March 1952, p. 68-73. 


This is an historical account of the development of these machines 
beginning with the period of GEORGE BooLe, the English scientist who first 
conceived the idea of symbolic logic as a system designed for the more 
efficient handling of problems. The paper discusses one of the devices of 
LULL, a 13th century Spanish mystic, who probably was the first to con- 
struct a machine to handle logical problems. The device consisted of a set of 
wheels having a common center. Around the edge of the disks were letters 
representing various ideas. When the wheels were rotated, various combina- 
tions of ideas would appear for further investigation. This primitive machine 
prompted CHARLES STANHOPE in the 18th century to invent the Stanhope 
demonstrator for solving syllogisms. Later WiLL1amM STANLEY JEVONS 
developed his “‘abecedarium,” a method of applying Boole’s approach. In 
the case of a syllogism all possible combinations of the three terms A, B and 
C are listed (the negative of these terms being represented by a, b, and c). 
In solving a problem all classes inconsistent with the premises are crossed 
out—leaving only the true combinations. These are then examined to deter- 
mine the relationship between the terms in an effort to weed out the false 
conclusions. Jevons went on to construct his more rapid “logical piano” on 
the same principle. 
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Since 1900 many devices for solving syllogisms have been invented. With 
the publication of an article by CLAUDE E. SHANNON of M.I.T. entitled, “A 
symbolic analysis of relay and switching circuits,”’ it was proved that circuits 
could be constructed in series and in parallel making use of the binary system 
to indicate the two truth values which correspond to basic Boolean relations. 
The Kalin-Burkhart Calculator, a small low-cost device, is the first in history 
to use electrical circuits for solving problems in symbolic logic. It is similar 
in principle to the Jevons piano but excells in that it can handle a greater 
number of problems far more intricate in complexity than those solved on 
Jevons’ machine. 

At present these machines have limited application; probably their 
chief value lies in the fields of business procedures and of electronic com- 
puters. In the latter field, especially, decisions involving considerations of 
logic are often necessary in setting up problems to be solved on the machine. 
It has even been suggested that with built-in probability values the machine 
could after considering all possibilities make important logical decisions. 
The question of whether or not such a machine could develop a creative 
imagination has been the subject of much fantastic speculation. 

Evita T. Norris 
NBSMDL 


1001. Joun D. GooDELL, ‘‘The foundations of computing machinery,” The 
Journal of Computing Systems, v. 1, June 1952, p. 1-13. TeEnnNy Lope, 
“The realization of a universal decision element,” The Journal of Com- 
puting Systems, v. 1, June 1952, p. 14-22. 


These articles are, according to their authors, the beginning of a series 
discussing the logical design of computing machinery in terms of Boolean 
algebra or, more generally, symbolic logic. They treat mainly the two- 
valued functions of one and of two two-valued variables, and particularly the 
recognition and realization of complete subsets of the two-valued functions. 
Functions of many arguments, variable functions and many-valued logics 
are very briefly discussed. 

R. D. ELBourN 
NBSEC 


1002. HARVARD UNIVERSITY, COMPUTATION LABORATORY, Annals, v. 27: 
Synthesis of Electronic Computing and Control Circuits. Cambridge, Mass., 
Harvard University Press, 1951, 278 p., 20 X 27 cm. Price $8.00. 


The bewildering complexity of electronic computing and control circuits 
and the profusion of means available for achieving the same functional 
performance make it extremely desirable to have 1) a concise and unique 
notation for writing the function of a circuit so that functionally identical 
requirements or circuits can readily be recognized, 2) a simple translation of 
this unique notation into a non-unique one corresponding more closely to the 
actual circuits so that one has a facile shorthand for manipulating a func- 
tional requirement into various possible realizations and 3) a method for 
finding a most economical realization of a required function. 

It is well known that Boolean algebra is useful for 1) and 2), but this book 
is the most extensive treatment of these applications which has appeared. 
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The method of “minimizing charts” introduced in chapter 5 is a significant 
contribution to the very formidable problem of most economical synthesis. 
At the present state of the theory, minimizing charts and tables included in 
the book provide a rather complete solution for circuits with one output and 
four or five input variables. For six or seven variables the charts become 
quite unwieldy. Some useful methods are given for circuits with multiple 
outputs and circuits containing signal delay or storage devices, but there 
remains the necessity for much arbitrary decision in their synthesis, and the 
economy problem is wide open. 

The second half of the book covers systems for binary coding decimal 
digits and the design of a variety of adders and multipliers mostly of the 
coded decimal type. 

The Boolean algebra of this book uses only the logical product X Y and 
negation, written X’ = 1 — X, because the arithmetic of these operations 
coincides with ordinary arithmetic. While this is a complete system, the 
reviewer would have included the logical sum, in the belief that the difficulty 
of remembering 1 + 1 = 1 is outweighed by having symmetrical notations 
for both types of rectifier circuits and by having the powerful duality prin- 
ciple to manipulate expressions; however this is merely personal preference. 

Although the efficient designing of electronic computing and control 
circuits has only begun to be systematized, anyone who wishes to use the 
best known techniques and possibly to contribute some new ones must 
certainly consult this volume. 


R. D. ELBourN 
NBSEC 


1003. OrFICE OF NAVAL RESEARCH, Digital Computer Newsletter, v. 4, July 
1952, 9 p. 


The contents are as follows: 


. Whirlwind I 
. The SEAC 
. Naval Proving Ground Calculators 
. Moore School Automatic Computer (MSAC) 
. The Logistics Computer 
. Aberdeen Proving Ground Computers 
The ENIAC 
The EDVAC 
The ORDVAC 
7. The Circle Computer 
8. Electronic Computer Corporation (ELECOM—100) 
9. The CADAC (CRD 102) 
0 
1 


10. The ACE Pilot Model, England 
11. The Ferranti Computer at Manchester University, England 


Data Processing and Conversion Equipment 


1. Digital to Analogue Converters 
2. Oscillograph Trace Reader 
3. Flying Typewriter 
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1004. INSTITUTE OF RADIO ENGINEERS, “Standards on Electronic Com- 
puters: Definitions of Terms, 1950.” Institute of Radio Engineers, Proc., 
v. 39, no. 3, 1951, p. 271-277. 


This is the first report of the IRE Subcommittee on Definitions of 
Electronic Computer Terms. R. SERRELL was chairman of thissubcommittee, 
J. W. ForreEsTER headed the Electronic Computers Committee, and J. G. 
BRAINERD was chairman of the Standards Committee. There is a total of 
approximately 135 terms for.which definitions are given. About 30% of these 
terms are arithmetical or logical, over one-half relate specifically to digital 
devices, and approximately 20 terms concern analog devices. (It might be of 
interest to note that a new IRE subcommittee is at present reviewing these 
definitions.) 

H. D. HuskEey 
Wayne University 
Detroit, Michigan 


1005. IBM Scientific Computation Forum, Proc. edited by H. R. J. Groscu, 
1948, 126 pages. 22 X 28.5 cm. 


This volume refers to one of a series of meetings arranged from time to 
time by the IBM Corp. (cf. the- following reviews). The contents are as 
follows: 


1. “Evaluation of higher order differences on the Type 602 Calculating 
Punch” by FRANK M. VERzUH 

2. “Differencing on the Type 405 Accounting Machine’ by GERTRUDE 
BLANCH 

3. “The use of optimum interval mathematical tables” by H. R. J .Groscu 
4. “Punched card techniques for the solution of similtaneous equations and 
other matrix operations” by WILLIAM D. BELL 

5. “Two numerical methods of integration using predetermined factors” by 
LELAND W. SPRINKLE 

6. ‘Integration of second order linear differential equations on the Type 
602 Calculating Punch’”’ by N. ARNE LINDBERGER 


7. “Integration of the differential equation —> = P-F(r) using the Type 


dr® 
601 Multiplying Punch” by PauL HERGET 
8. “‘Some elementary machine problems in the sampling work of the census” 
by A. Ross ECKLER 
9. “IBM applications in industrial statistics” by CUTHBERT C. HurD 
10. “Some engineering applications of IBM equipment at the General 
Electric Company” by FRANK J. MAGINNISS 
11. “Planning engineering calculations for IBM equipment” by BEN 
FERBER 
12. ‘A survey of the IBM project at Beech Aircraft Corporation” by JoHN 
KINTAS 
13. “‘Aerodynamic lattice calculations using punched cards” by HANS KRAFT 
14. “Dynamics of elliptical galaxies’’ by Jack BELZER, GEORGE GAMOW 
and GEOFFREY KELLER 
15. ‘Application of punched cards in physical chemistry” by GILBERT W. 
KING 
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16. “Application of punched card methods to the computation of thermo- 
dynamic properties of gases from spectra” by Lypia G. Saveporr, Jack 
Belzer and H. L. JoHNsTon 

17. “Calculation of the equilibrium composition of systems of many con- 
stituents” by StuaRT R. BRINKLEY, JR. and RoBERT W. SMITH, JR. 

18. ‘Punched card calculating and printing methods in the Nautical 
Almanac Office’ by FREDERICK H. HOLLANDER 

19. “Programming and using the Type 603-405 Combination Machine in 
the solution of differential equations’ by GrorGE S. FENN 

20. “Applications of punched card equipment at the Naval Proving Ground” 
by Cuiinton C. BRAMBLE 

21. ‘‘Use of the IBM relay calculators for technical calculations at Aberdeen 
Proving Ground” by JosepH H. LEvIN 

22. ‘Simultaneous linear equations’”’ by FRANCIS J. MURRAY 

23. ‘Computation of shock wave refraction on the selective sequence 
electronic calculator” by HARRY POLACHEK 

24. “Computation of statistical fields for atoms and ions’”’ by L. H. Tomas 


At the meeting reported here, most of the papers described details of the 
use of IBM machines for various problems in numerical mathematics, which 
are adequately characterized by their titles. A few of the papers present 
surveys of the computing work done in certain organizations on IBM 
machines. The paper of Fenn introduces the 603-405 combination (cf. the 
following review). Murray's paper refers to his special-purpose machine for 
solving linear equations. Grosch’s article is of more special interest to table- 
makers. 


NBSCL 


1006. IBM Seminar on Scientific Computation, Proc. edited by CUTHBERT 
C. Hurp, November 1949, 109 pages. 22 X 28.5 cm. 


The contents of this volume are as follows: 


“The dynamics of nuclear fission” by Davin L. HILL 

. “Monte Carlo calculations” by W1rLL1am W. WoopsurRY 

“Modification of the Monte Carlo method” by HERMAN KAHN 
“Analyzing exponential decay curves’’ by ALSTON S. HOUSEHOLDER 
“On the distribution of Kolmogorov’s statistic for finite sample size’ by 
. W. BiRNBAUM 

. ‘The IBM Card-Programmed Electronic Calculator’? by CutHBert C. 
Hurp 

7. “Stochastic methods in quantum mechanics’ by GiLBERT W. KING 

8. “Calculation of resonance energies” by GEORGE E. KIMBALL 

9. “Cam design calculations on the Card-Programmed Electronic Cal- 
culator” by E. A. BARBER 

10. “Calculation of the equilibrium composition of homogeneous multi- 
component systems’”’ by Stuart R. BRINKLEY, JR. and Rosert W. 


AN 


R. 
11. “Eigenvalue problems related to the Laplace operator” by DonaLp 
A. FLANDERS and GEorGE H. SHORTLEY 
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12. “Numerical solution of partial differential equations of parabolic type” 
by L. H. THomas 

13. ‘Solutions of the wave equation” by Paut HERGET 

14. “Sampling methods applied to differential and difference equations” by 
Joun H. Curtiss 


The conference about which this volume reports was on a higher scientific 
level than other similar conferences arranged by the IBM Company and 
requires more detailed review. Not all the papers in this volume are primarily 
of interest to the numerical analyst. For example, Kahn's paper is an 
exposition of Monte Carlo methods as applied to definite integrals and inte- 
gral equations, with emphasis on importance sampling, but without reference 
to machine methods or to any of the deeper questions of numerical analysis. 
The main purpose of Birnbaum’s paper is “‘to point out that in view of the 
development of high-speed sequence computing equipment, the tabulation of 
the exact probability distribution of Kolmogorov’s statistic for finite N has 
become practically feasible, and to propose that such a tabulation should be 
carried out.” The article by G. W. King is principally concerned with the 
formulation of certain quantum mechanical problems in terms of finite- 
difference equations, and thence in terms of random walks; some particular 
numerical results are given as an indication of the rate of convergence of the 
Monte Carlo method when used for the solution of such problems. Similarly, 
Kimball’s paper is principally interesting for the mathematical formulation 
of a physical problem rather than for the method of solution. 

Other papers contain more substantial contributions to our understand- 
ing of machine methods of computation. The article by Hill outlines the 
theory of the liquid-drop model of.nuclear fission and gives details of the 
numerical techniques used on the SSEC for certain prototype computations 
in this field. The paper by Woodbury describes shielding computations by a 
Monte Carlo method, run on the ‘405-603 Combination,” the famous ma- 
chine improvised by Woodbury and Toben by combining two standard 
IBM machines. This combination was the predecessor of the IBM Card- 
Programmed Calculator. The approximation of a given function by a linear 
combination of exponentials is the subject of Householder’s paper. Several 
methods—graphical, least squares, iterations based on Poisson-distributed 
errors—are compared. C. C. Hurd gives a functional description of the 
Card-Programmed Calculator. This is the so-called ‘“‘Model I’’ which was in 
process of construction at the time of the presentation of this paper. Barber’s 
paper describes the evaluation of certain polynomials. Brinkley and Smith 
report on progress in their long-standing problem which consists mathemati- 
cally of the solution of numerous systems of nonlinear algebraic equations. 
They describe a special method applicable in certain degenerate cases and 
mention briefly the use of the Card-Programmed Calculator. The paper of 
Flanders and Shortley is a discussion of several iterative methods for the 
computation of wave functions; particularly methods employing matrix 
operators which are obtainable by forming polynomials of simpler operators. 
Herget gives a sketchy description of the numerical integration of a system of 
second order ordinary differential equations on an IBM 602-A calculator. 

The most important papers in this volume, from the standpoint of numer- 
ical analysis, are those of L. H. Thomas and J. H. Curtiss. The former is an 
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exhaustive treatment of questions of approximating parabolic differential 
equations by finite-difference methods, the errors in the solution caused by 
such approximation, stability and rate of convergence. While the discussion 
is purposely nonrigorous, it is fundamental and complete. An unpublished 
but frequently quoted method due to J. von NEUMANN is used. 

The paper by J. H. Curtiss, which occupies almost one-fourth of the 
volume, is the only one written with the care and rigor customary in mathe- 
matical publications, while the rest of the volume is in the nature of an 
exchange of information on techniques, with emphasis on timeliness, 
effectiveness in getting results, and “know-how” rather than on scientific 
requirements. The paper begins with a historical survey of the field, which is 
shown to be much older than its currently fashionable name of “‘Monte 
Carlo Method.” There follows an exposition which brings together results, 
some of which had previously been widely scattered and relatively inacces- 
sible in the literature on mathematical physics and theory of probabilities. 
Significant new results are presented on the standard error of the solution, 
the number of samples (random walks) necessary to achieve a given accuracy, 
the mean length of a walk and thus the expected amount of computing effort 
necessary for a given problem, and on importance sampling. 

The editorial work is excellent throughout the volume, except for the 
discussions at the end of the papers, where misprints like the following occur: 
p. 27, 1.6, for “‘Prof. Katz and Dr. Dunning”’ read “‘[M.] Kac and [M. D.] 
DonSsKER”’; p. 78, 1. 22, for “‘Dr. Foster’ read “‘[G. ForsyTHE.”’ 

F. L. ALT 
NBSCL 


1007. IBM Industrial Computation Seminar, edited by CUTHBERT C. Hurp, 
September 1950, 103 p. 22 XK 28.5 cm. 


Participants in this seminar discussed fundamental computational meth- 
ods used by research engineers and scientists in a wide variety of research 
problems. Particular attention was drawn to computational techniques 
developed in the fields of chemistry and petroleum. The contents of this 
volume are as follows: 


1. “The role of the punched card in scientific computation” by WALLACE J. 
ECKERT 

2. “Machine calculation of the plate-by-plate composition of a multicom- 
ponent distillation column” by ASCHER OPLER and RoBeErt C. HEITz 

3. “Continuous distillation design calculations with the IBM Card-Pro- 
grammed Electronic Calculator’ by ARTHUR ROSE, THEODORE J. WILLIAMS 
and S. Dye, III 

4. “Application of automatic computing methods to infrared spectroscopy” 
by GILBERT W. KING 

5. “Correlation and regression analysis’”’ by E. L. WELKER 

6. ‘‘Pile-driving impact’’ by Epwarp A. SMITH 

7. “Punched card mathematical tables on standard IBM equipment’’ by 
ELEANOR KRAWITZ 

8. “The solution of simultaneous linear equations using the IBM Card- 
Programmed Electronic Calculator” by Justus CHANCELLOR, JOHN W. 
SHELDON and G. Liston TATUM 


17 
M 
18 
ar 
G 


B@eanggs 


Z 


9. 

TI 
ie 

Ay 

10 

Pr 

11 

13 

an 

14 

EI 

2 R 

16 
te 
“we 

= 


rd- 


AUTOMATIC COMPUTING MACHINERY 55 


9. ‘Two applications of the IBM Card-Programmed Electronic Calculator: 
The Gauss-Seidel method of solution of simultaneous linear equations; 
Approximating the roots of a polynomial equation” by I. C. Liccett 
10. ‘Matrix by vector multiplication on the IBM Type 602-A Calculating 
Punch” by ELEANOR KRAWITZ 
11. “Numerical solution of two simultaneous second-order differential 
equations” by WALTER H. JOHNSON 
12. “Numerical evaluation of integrals of the form S8f(x)g(x)dx” by Joun 
W. SHELDON 
13. “‘The use of orthogonal polynomials in curve fitting and regression 
analysis” by JACK SHERMAN 
14. “‘General-purpose ten-digit arithmetic on the IBM Card-Programmed 
Electronic Calculator” by Stuart R. BRINKLEY, Jr., G. L. WAGNER and 
R. W. JR. 
15. “‘Remarks on distillation calculations’’ by Joun W. DoNNELL 
16. Some applications of the Monte Carlo method: 
“Matrix inversion on the IBM Accounting Machine” by Ascher 
Opler 
“Remarks on finding roots of, and inverting, a matrix” by GILBERT 
W. KInG 
“Remarks on the Monte Carlo method” by Cutapert C. Hurp 
17. “Plotting punched card data using the IBM Type 405 Accounting 
Machine”’ by Paut T. Nims 
18. “A method for evaulating determinants and inverting matrices with 
arbitrary polynomial elements by IBM punched card methods” by L. E. 
Grosu, Jr. and E. Upsin 
The article by Chancellor, Sheldon and Tatum describes a method for 
solving simultaneous linear algebraic equations up to and including order 
twenty-one using the Crout method. With this procedure, a 20 X 20 matrix 
can be inverted in two hours. The prescribed method assumes a model I 
CPC equipped with one (941) auxiliary storage. The article by Brinkley, 
Wagner and Smith, Jr., describes a general purpose ten-digit arithmetic set- 
up for the model I CPC. In addition to the elementary arithmetic operations, 
square root and certain combinations, one can automatically select alter- 
native computational routines and also select input data. 
R. K. ANDERSON 
NBSCL 


1008. IBM, Industrial Computation Seminar, Proc. edited by Cuthbert C. 
Hurp, August 1951, 148 p. 28.5 K 22 cm. 


The participants of this seminar were research engineers and scientists 
representing computing facilities which employ IBM card-programmed 
electronic calculators. The contents of this volume are as follows: 


1. “Application of the IBM Card-Programmed Electronic Calculator to 
engineering procedures at the Glenn L. Martin Company” by W. B. Kocu 
2. “Reduction of six-component wind tunnel data using the IBM Card- 
Programmed Electronic Calculator, Model II’ by M. L. LEssER 

3. “IBM Card-Programmed Electronic Calculator operations using a Type 
402-417BB and 604-2” by H. E. Tituitt, M. Kenyon and B. OLDFIELD 
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4. “The Combomat” by J. D. MADDEN 

5. “The IBM Type 604 Electronic Calculating Punch as a miniature elec- 
tronic calculator’ by P. T. 

6. “General-purpose floating point control panels for the IBM Card- 
Programmed Electronic Calculator” by N. A. Patton, K. BERGER and L. R. 
TURNER 

7. “Catapult takeoff analysis” by J. R. Lowe 

8. “Computation of loan amortization schedules on the IBM Card-Pro- 
grammed Electronic Calculator” by C. H. GuSHEE 

9. “Techniques for handling graphical data in the IBM Card-Programmed 
Electronic Calculator” by W. D. BELL 

10. “Calculation of the flow properties in an arbitrary two-dimensional 
cascade” by J. T. HORNER 

11. “Automatic calculation of the roots of complex polynomial equations 
using the IBM Card-Programmed Electronic Calculator” by J. GALLISHAW, 
Jr. 

12. “A recursion relation for computing least square polynomials over 
moving-arcs” by GEORGE R. TRIMBLE, JR. 

13. ‘Numerical solution of second-order non-linear simultaneous differential 
equations” by Henry S. WOLANSKI 

14. “Matrix inversion and solution of simultaneous linear algebraic equa- 
tions with the IBM Type 604 Electronic Calculating Punch’’ by GrorGe W. 
PETRIE, III 

15. “The determination of eigenvectors and eigenvalues of symmetric 
matrices” by EVERETT C. YOWELL 

16. “An application of the IBM Card-Programmed Electronic Calculator to 
analysis of airplane maneuvering horizontal tail loads’’ by LoGAN T. WATER- 
MAN 

17. “Fifth-order aberration in an optical system” by Ruta K. ANDERSON 
18. “Theory of plastic vibrations of helicopter fuselages” by PETER F. 
LEONE 

19. ‘Machine procedure for computation of elastic vibrations of helicopter 
fuselages’”’ by WILLIAM P. HEISING 


Nims describes how a 604 Electronic Calculating Punch may be used as a 
small sequence-controller calculator. The 604 can perform a different arith- 
metic operation on each card. In addition it is possible for the machine to 
shift, look up entries in a table of functions, calculate square roots and stop 
when certain check conditions do not hold. 

Patton, Berger and Turner have designed general-purpose floating point 
control panels for.the CPC which give the coder a great deal of flexibility. 
In addition to the basic arithmetic operations, one has three independent 
program fields and four numerical fields (per card), automatic extraction of 
integral roots from second to seventh, eight kinds of conditional operations 
and machine stop with automatic list cycles on all improper operations 
(division by zero, overflow). 

Gallishaw’s technique for finding complex roots of a polynomial is an 
iterative procedure employing synthetic division and Newton’s method. 
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Using standard trial values, it takes from twenty to twenty-five minutes to 
find all the roots of a complex 7th degree equation. If any of the roots are 
approximately known in advance, they can be used instead of the standard 
trial values, thus cutting down the solution time. The process uses eight- 
digit arithmetic throughout. 

Petrie has devised a method for matrix inversion and solution of simul- 
taneous linear equations using a 604 Electronic Calculating Punch, a tabu- 
lator and two reproducers. 

The procedure consists of a continuous flow of cards through these four 
machines, the cycle being repeated N times to invert an Nth order matrix 
with approximately N? cards each time. The size of the matrix is unlimited, 
and a check operation is included. 

R. K. ANDERSON 
NBSCL 


1009. IBM Technical News Letter No. 3, Applied Science Department, 
December 1951, 106 pages, 22 X 28 cm. 


The contents are as follows: 


1. ‘Mass spectrometer calculations on the IBM Type 602-A Calculating 
Punch” by W. H. Kine, Jr. and WILLIAM PRIESTLY, JR. 

2. “Computations of inverse matrices by means of IBM machines” by 
JAacK SHERMAN 

3. “Linear equations and matrix inversion” by Eric V. HANKAM 

4. “Matrix and vector algebra” by Eric V. Hankam 

5. “Notes on the IBM Type 604 Electronic Calculating Punch and Type 
602-A Calculating Punch” by Eric V. Hankam 

6. “Stop controls for the IBM Type 604 Electronic Calculating Punch,” 
by BrusE MONCREIFF 

7. “An eight-digit general-purpose control panel” by WILLIAM P. HEIsSING 
8. “Interpolation on the IBM Card-Programmed Electronic Calculator” by 
Stuart R. BRINKLEY, JR. and G. L. WAGNER. 


The article by King and Priestly describes in detail the control panels 
used in performing a matrix multiplication of the 20th order and in sub- 
sequent normalization steps. The general-purpose board described by Heising 
(for use with either the IBM Type 604 Electronic Calculator or the IBM 
Card-Programmed Calculator) calculates, in addition to the basic arithmetic 
operations, VA, exp (+ A), log (1 + A), sin A, cos A, sinh A, cosh A, arcsin 
A, arctan A, arcsinh A, and arctan A. 

R. K. ANDERSON 
NBSCL 


NEws 


Association for Computing Machinery.—The University of Toronto Computation 
Centre acted as host to the Association for the fall meeting held September 8 and 9, 1952. 
The University had on display for its guests the three computers located at the Computation 
Centre, namely, FERUT, a newly installed large-scale computer built by Ferranti Limited 
which was undergoing tests at the time; UTEC, the model computer built at the University; 
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and BERTIE, a noughts and crosses machine built by Rogers Majestic Electronic Labora- 
tories. The program for the meeting was as follows: 


Sept. 8, 1952 10:00 a.m. to 12 noon 


General Session S. B. Witt1aMs, President, ACM, Chairman 
Welcoming address S. SmitH, President, Univ. of Toronto 
Pure and Applied programming M. V. Wi1LKEs, Director, University Mathe- 

matical Laboratory, Cambridge, England 
MANIAC N. Metropouis, E. F. Krein, W. Orve- 


DAHL, J. R. RicHarpson, H. B. Demuts, 
J. B. Jackson, Los Alamos Scientific 
Laboratory 


Chebyshev polynomials in the solution of | C. Lanczos, NBSINA 
large scale linear systems 


Sept. 8, 1952, 2:00 p.m. to 5:00 p.m. 


Programming C. C. Hurp, IBM, New York, Chairman 
Interpretative sub-routines J. M. Bennett, D. G. Prinz, M. L. Woops, 
Ferranti, Ltd., Manchester, England 
Machine aids to coding E. J. Isaac, NRL 
Computer aids to code checking I. C. Dresw, NBSCL 


Input scaling and output scaling for a E. F. Copp, H. L. Herricx, IBM 
binary calculator with decimal output 
Compiling routines R. K. Rmpcway, Eckert-Mauchly Division 
of Remington RAND, Inc., Philadelphia 
Function table techniques in programming R. F. SHaw, Electronic Computer Corp., 
New York 


Sept. 8, 1952, 2:00 p.m. to 5:00 p.m. 


Computing machines K. F. Tupper, Dean, Faculty of Applied 
Science and Engineering, Univ. of 

Toronto, Chairman 
The logical design of the Oak Ridge digital C. L. Perry, Oak Ridge National Labora- 


computer tory 

The Oak Ridge automatic computer J. C. Cau, Argonne National Laboratory 

Designing a low cost general-purpose W.E. Dossins, Computer Research Corp., 
computer Hawthorne, California 

The University of Toronto model elec- R. F. JonNnston, Computation Centre, 
tronic computer Univ. of Toronto 

A description of the electronic computer G. Estrin, Institute for Advanced Study, 
at the Institute for Advanced Study Princeton 

Banquet, Sept. 8, 1952, 7:30 p.m. 
Guest Speaker 
The St. Lawrence Seaway Project R. H. Saunpers, Chairman, Hydro-Elec- 


tric Power Commission of Ontario 


Sept. 9, 1952, 9:00 a.m. to 12:00 noon 


Machine calculations Mina Rees, ONR, Chairman 
Errors in iterative solutions of linear sys- A. S. HousEHOLDER, Oak Ridge National 
tems Laboratory 
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The numerical solution of a partial differ- 
ential equation on the IBM Type 701 
Electronic Data Processing Machines 

A numerical solution of the helium wave 
equation with the SEAC 

Matrix inversion by partitioning 


The use of subroutines on SEAC for nu- 
merical integrations of differential equa- 
tions and for Gaussian quadrature 

Use of continued fractions in high-speed 

computing 


D. W. Lapp, J. W. SHetpon, Applied 
Science Dept., IBM, New York 


J. H. Wecstern, NBSCL 
M. Lorxin, R. RemMaGe, BRL, Aberdeen 


Proving Ground 
P. Rasrnowitz, NBSCL 


D. Te1curoEw, NBSINA 


Sept. 9, 1952, 9:00 a.m. to 12:00 noon 


Storage devices 
The testing of cathode ray tubes for use in 
the Williams type storage system 
Williams tubes selection program 


Improvement of Williams memory relia- 
bility 

Ferro electric materials as storage ele- 
ments for digital computers and switch- 
ing systems 

An experimental rapid access memory 
using diodes and capacitors 


J. A. Raycuman, RCA, Chairman 

A. Rostinson, Ferranti, Ltd., Manchester 
England 

J. C. Cuu, R. J. Kier, Argonne National 
Laboratory 

R. Scnumann, Argonne National Labora- 
tory 

J. R. AnpErRson, BTL, New York 


A. W. Hott, NBSCML 


Sept. 9, 1952, 2:00 p.m. to 5:00 p.m. 


Logical design 


Symbolic synthesis of digital computers 
Logical or non-mathematical programmes 


A simplified universal Turing machine 
Simple learning by a digital computer 


An analysis by arithmetical methods of a 
calculating network with feedback 


W. H. Watson, Head Dept. of Physics, 
Univ. of Toronto, Chairman 

I. S. Reep, MIT 

C. S. StracHey, National Research Devel- 
opment Corp., London, England 

E. F. Moore, BTL, New Jersey 

A. G. OETTINGER, The Computation Labor- 
atory, Harvard Univ. 

L. C. Rossrns, Burroughs Adding Ma- 
chine Company, Philadelphia 


Sept. 9, 1952, 2:00 p.m. to 5:00 p.m. 


Engineering 
A high-speed magnetic-core output printer 


Development of computer components 
and systems 


An electronic analogue machine for com- 
puting the roots of algebraic equations 
of degrees through the eighth 
Analogue calculation of polynomials 
and their zeros 


Operating efficiencies and characteristics 
of the computing machines at Aberdeen 
Proving Ground 

Installation of a large electronic computer 


B. V. Bowben, Ferranti, Ltd., Manchester, 
England 

B. M. Gorpon, R. N. Nicora, Laboratory 
for Electronics, Inc., Boston 

W. S. Etxiott, Research Laboratories of 
Elliott Brothers (London), Ltd., Bore- 
hamwood, England 

L. LorGREn, Research Institute of National 
Defense, Radio Department, Stockholm, 
Sweden 

M. G. SCHERBERG, J. F. Rrorpan, Wright 
Air Development Center, Wright-Patter- 
son Air Force Base, Ohio 

H. Spence, Ballistic Research Laboratories, 
Aberdeen Proving Ground 


L. R. Jounson, Hdgs., U. S. Air Force 
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NBSNAML.—On May 15, 1952, at Washington, D. C., the National Applied Mathe- 
matics Laboratories sponsored a meeting of Mathematical Tables in the Light of Electronic 
Computers. The program for the meeting was as follows: 


Work in Progress J. Topp, NBSCL, Chairman 
E. T. Goopwin, Royal Society 
Mathematical Tables Committee 
J. P. Wonc, RAND Corporation 
M. NBSCL 
G. BLancn, NBSINA 


(There was also a written communication from M. PiconE, Instituto Nazionale per le 
Applicazione del Calcolo, Rome.) 


Current Needs ~ F. L. Att, NBSCL, Chairman 

Physics G. Breit, Sloane Physics Dept., Yale Univ. 

Aerodynamics C. W. Jones, Dept. of Math., Liverpool 
Univ. 

Military applications L. S. DEDERICK, Ballistic Research Labora- 
tories, Aberdeen Proving Ground 

Statistics C. and D. TEICHROEW, 
NBSSEL 


Is the preparation of a single volume 7- R. C. ARCHIBALD, Brown Univ. 
place table of J,(x) of integral order up 
to 120, at interval .01, desirable? 

Are tables of Bessel functions of thesecond R. C. ARcHIBALD, Brown Univ. 
kind and of orders higher than 20 de- 
sirable? 


Long Term Policy P. M. Morse, MIT, Chairman 
Round table discussion of which the principal speakers were: 
C. C. Hurp, IBM Corporation 
D. H. LenmMer, NBSINA 
E. T. Goopwin, Royal Society 
Mathematical Tables committee 
J. H. Curtiss, NBSNAML 


Editorial Matters D. H. Lesmer, NBSINA, Chairman 
Index of Statistical Tables D. L. Wattace, Dept. of Mathematics, 
Princeton Univ. 
Index of Mathematical Tables C. W. Jones, Dept of Mathematics, 
Liverpool, University 
Editing of Mathematical Tables G. M. CLemMEeNcE, Naval Obsrevatory, 


Washington, D. C. and D. H. SapLer, 
His Majesty’s Nautical Almanac Office 
The preparation of a bibliography of cal- R.C. ArcuIBALD, Brown Univ. 
culating machines developed from lists 
of Aiken, Travis and an Australian 
physicist, from Patent Office records, 
and from other sources. 
An organization to assemble for MTAC _iR.. C. ARcHIBALD, Brown Univ. 
tabular material news from the very 
numerous research centers of this and 
other countries 
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OTHER AIDS TO COMPUTATION 


ANALOGUE CALCULATION OF POLYNOMIAL AND 
TRIGONOMETRIC EXPANSIONS 


Introduction: The use of polynomials and trigonometric expansions in 
science and engineering is well known, and hence the importance of machine 
application for their rapid evaluation does not need strong defense. In fact, 
the multiplicity of special purpose polynomial computers that have been de- 
signed and built points to the great desire for machine treatment of these 
functions. Although the method described in this paper is applicable to dig- 
ital as well as analogue machines, the original intent of the authors was 
slanted towards the use of the analogue machine, and hence the presentation 
has an analogue machine background. The significance and importance of 
the new method is its simplicity and its use of standard computing machine 
equipment. 

The Vector Multiplier: The key process which leads to the new method 
is the automatic multiplication of an arbitrary pair of vectors. The present 
paper will be limited to the treatment of the complex vector (a + 1b). The 
product of two such vectors, 2; = x; + iy; and 2: = x2 + %y2, may be 
written 


(1) = 2182 = — + + X21), 
and the operation performed by machine as shown in Figure 1. 


J. 


sign-changer sumers 


‘ 


2 


23" (x,25- ¥,¥2) 1 


Fic. 1. The Vector Multiplier. 


In terms of operation on either the Reeves REAC or the comparable 
Goodyear Electronic Machine, x: and y; in the figure represent input 
voltages to a pair of ganged potentiometers or a set of four simple potentiom- 
eters (hereinafter referred to as pots) set respectively at x2 and ye. The 
outputs of these pots are, except for sign, the four terms on the right of 
equation (1). The term yxy is put through a sign changer, and then the terms 
are put into a pair of summers to obtain outputs which are the negatives of 
the real and imaginary parts of z3. The amplifiers used on the machines re- 
ferred to automatically change the sign of their summed inputs and hence 
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we find the complex vector multiplier described above giving an output 
vector which is the negative of the product 23. This is in no way an impedi- 
ment in the operation since many operations will call for sign changes, and 
hence, on the average the same sign change equipment would be needed even 
if the vector product operation produced the product without sign change. 
In what follows we shall call the input z; the multiplicand and the pot setting 
the multiplier. 

Polynomials in the complex variable z with real coefficients: It is now 
fairly obvious how successive application of the vector multiplication will 
produce the successive powers 2’, 2’, 2, ---, 2". Figure 2 shows a part of the 
circuit diagram for a fifth degree polynomial. Unit voltages enter at the left 
and the successive output voltages representing respectively 2, — 2’, 2°, 
— 2 and 2° are shown by the vertical output lines. If ganged pots are avail- 


§ 


Fic. 2. Circuit which produces the real and imaginary parts of z, — 2*, 2*, — 2‘ and 2°. 
= x + iy) 


able, then the x pots and the y pots are set simultaneously by two or more 
operations depending on the multiplicity of pots per control. Of course the 
servo ganged pots usually provided with analogue calculators may be used 
for the ganged pot operation. Figure 3 shows the servo-pot diagram equiva- 
lent of Figure 2, with additional sign changers to provide the powers of z 
with both signs. 

The outputs R(z), — R(z*), R(z*), — R(z*) and R(z‘), i.e., the real parts 
of the successive powers indicated, are respectively put through a set of pots 
which multiply them by the polynomial coefficients a1, a2, a3, ---, as. The 
outputs of the coefficient pots, including a pot representing the polynomial 
constant, with appropriate signs, then go into a summing amplifier whose 
output defines the real part of the polynomial or its negative depending upon 
the computing requirements. A comparable process is used to calculate the 
imaginary part of the polynomial. 

Polynomials in z with complex coefficients: In this instance it is only a 
matter of making vector multiplications of the coefficients a, into the 
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corresponding vector terms R(z”) + iI(z™) in place of using simple pot 
multiplication. The summing amplifiers used in the latter multiplications 
may of course be used for all or part of such multiplications depending upon 
the amplifiers input facilities. The real and imaginary parts of the polyno- 
mial may be projected on the screen of the cathode ray tube so that a visual 
display of the calculated points becomes available. Through this display one 
discovers rather quickly characteristics of the polynomial under study such 
as dominant terms and zeros. 


+ 


Fic. 3. Circuit in which the potentiometers of Fig. 2 have been replaced 
by servo multipliers. 


Trigonometric forms in one variable: Consider forms reducible to 


(2) T(0) = ao + a1 cos 6 + a2 cos 26 + --- +a, cosn@ 
+ sin 6 + b.sin 20 + --- +5, sin 


in which aj, a2, Gn, 6, bs, are constants and m, are integers with 
m > n. The form (2) is the real part of the following polynomial : 


(3) ao + (a + (a2 + + (a, 
+ + Gn — 


in which Ga41, @n42, ***, @m are all zero. To evaluate the trigonometric form 
for any specific range of 6 one has but to evaluate the polynomial (3) over an 
x, y range corresponding to cos @, sin @. 

The resolver servos available on the machines referred to would facilitate 
such treatment of the independent variables. In any case, the use of the sine 
and cosine columns in the trig table would make it relatively simple to range 
the variable appropriately. The treatment of mixed expansions of terms of 
the type r* cos 28, r* sin ”@ is obvious. 
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Polynomials in Several Variables: We consider first a polynomial 
P(21, 2) in two variables (z:, 22). In this case the successive powers in 2; and 
Z2 are respectively formed by the method described above and the required 
cross products such as 2;", 2" are formed by the required vector multiplica- 
tion. The vector cross product terms are then multiplied by the appropriate 
scalar or vector coefficient and the results fed into summing amplifiers to 
form the real and imaginary parts of P(z;, 22). Treatment of the more general 
case is obvious. In the case of polymonials in several variables this machine 
calculation becomes very useful for obtaining exploratory information about 
the characteristics of such functions. 

Polynomial transformations, approximations, zeros: When a polynomial 
has been coded on a machine, it becomes possible to obtain information by 
either manipulating the coefficient pots or the z pots or both. By using ca- 
thode ray screens and function tables one manipulates the coefficients to 
determine a transformation to take a contour C in the z plane into a contour 
C’ in a P plane or one finds a polynomial approximation to a given contour 
C’ in a P plane. Again one manipulates the x, y pots and determines a set of 
points in the z plane corresponding to a point p in a P plane, in other words 
one solves the equation p = P(z) for a given value of p. 

A systematic process for treatment of the latter problem is known in the 
literature’? and will here be outlined briefly for the sake of completeness. 
The brief statement will enable most readers to apply the process without 
referring to the literature. We have the theorem that a closed contour C 
in the z plane is mapped by a polynomial inte a closed contour C’ enclosing 
a given point » if and only if C contains at least one of the points 2, satis- 
fying the equation p = P(z,). One finds a rectangle C containing a point z, 
and then gradually contracts it by systematic trial and error until z, is 
located with the required accuracy or within the accuracy of the machine. 
Further accuracy may be obtained by the complementary use of a hand digi- 
tal machine by such method as found in MILNE.’ 

It must be pointed out again that although the machine coding has been 
described in terms of ganged pots the coding may also be done with simple 
pots when ganged pots are not available. In fact the example cited below 
was worked out with and without ganged pots. With simple (x, y) pots 
all of which must be reset for each change in z, the treatment of the problem 
is somewhat longer. It was observed that some of the simple x, y pots 
dominate the changes in the polynomial and with this experience it did not 
take long to drive the polynomial values in the required direction. 

Example: The polynomial equation 


P(z) = — — 32° +27 +122+5=0 


was solved on both Reeves and Goodyear Analogue Machines. One root was 
located in the rectangle whose vertices were (0, 0), (2, 0), (2, 27) and 
(0, 24). The rectangle was systematically contracted until the root was 
located with an optimum of accuracy. The polynomial was then reflected in 
the origin and a pair of roots located in the same rectangle. Since one of 
these was clearly a real root it was located with an optimum of accuracy 
by an obvious systematic search on the real segment of the rectangle. The 
first contraction of this rectangle which left out the axis of reals gave a con- 
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tour in the P plane which simply inclosed the origin. It was then a matter of 
further contractions to obtain the best estimate of the root. The table below 
gives a comparison of the roots as calculated by both analogue and digital 
methods. 


Analogue Method — .454 1.94 + 9601 — 1.21 + .954 
Digital Method — 4517 1.935 + .9466i — 1.216 + .95954 


Max G. SCHERBERG 
Joun F. R1oRDAN 
Wright-Patterson Air Force Base 
Ohio 
This paper was written at the Flight Research Laboratory, Wright Air Development 
Center, Wright-Patterson Air Force Base, Ohio. 
1L. BAvER & S. Firer, “The solution of polynomial equations on the REAC,” Sympo- 
sium I on REAC Techniques. U. S. Navy Special Devices Center and Reeves Instrument 
New York, 1951, p. 31-36. 
od TrrcHMaRsi, Theory of Functions. Oxford University Press, 1939, Second 
tion, 
Mune, Numerical Calculus. Princeton University Press, 1949, p. 53-57. 


BIBLIOGRAPHY Z 


1010. J. Euter & R. Lupwic, “Zwei Nomogramme zum Gebrauch bei 
Messungen mit optischen Pyrometern,” Zeitschrift f. angew Physik, 
v. 2, 1950, p. 362-373. 


The nomograms involve two-dimensional plotting for the inputs. 
F. J. M. 


1011. H. Gtusrecnt, “Elektrisches Rechengerat fiir Gleichungen héheren 
Grades,” Zeitschrift f. angew. Physik, v. 2, 1950, p. 1-8. 


The apparatus is set up to represent both the z and w plane and the roots 
of the polynomial are located by exploration. The null point of the poly- 
nomial is obtained by means of a special cathode ray tube arranged so that 
when both the real and imaginary parts of w are zero an eiectron beam goes 
through a narrow hole and signals the presence of a zero. As far as practical, 
quantities are represented by voltages having complex significance, thus 
multiplication is by means of a mixing tube and the resolution into real and 
imaginary parts is accomplished by a passive network associated with the 
ordinary type of cathode ray tube upon which z and w are represented. 
Normally a circle is mapped in the z plane and its w image also is given. 
The amplitude of this circle can be either varied by hand or swept auto- 
matically. The article discusses also two other apparatuses of this type due 
to TISCHNER and Rascu. 

F. J. M. 


1012. E. H. GAMBLE & B. W. HatTEn, “Design and analysis of a conserva- 
tive dynamic load simulator,” Jn. Appl. Phys., v. 22, 1951, p. 1250-1257. 


The design of a load simulator for the dynamic testing of actual control- 
system components such as aircraft autopilots is described. The use of an 
electrical computer for determining the load and a hydraulic servomecha- 
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nism for actually applying the load to the test unit results in greater flexi- 
bility than a load simulator in which the loading is determined mechanically. 
A feature of this simulator which should prove of great importance in the 
testing of large servomechanisms is the fact that virtually all the energy 
absorbed from the unit being tested is returned to the electrical lines and 
need not be dissipated as heat. 

RIcHArRD C. Booton, Jr. 
Massachusetts Institute of Technology 
Cambridge 39, Massachusetts 


1013. G. A. Korn & T. M. Korn, Electronic analog computers (D-c analog 
computers): McGraw-Hill Book Co., New York, 1952, xv + 378 p., 
15 X 23 cm. Price $7.00 


This book is concerned with d.c. analog equipment of the type manu- 
factured by Reeves, Goodyear and Boeing and to a certain extent with that 
of Philbrick. There is also a clear intention to indicate how less expensive 
versions of these machines may be set up for special purposes. Two special 
installations, the Curtiss Wright Analog Computer and the RAND analog 
facility, are discussed. 

The first three chapters discuss the principles of these devices, practical 
set up procedure and a variety of specific applications. The fourth chapter is 
concerned with the theory of linear computing networks and specifically the 
necessity for using amplifiers and the effect of this use on the solution. The 
fifth chapter discusses the design of d.c. amplifiers for computer applications. 
Various simple types of such amplifiers are discussed as well as the modern 
drift stabilized amplifiers as developed by Wi1LLiaMs, TARPLEY & CLARK, 
and the GOLDBERG amplifier. The titles of Chapters 6, 7 and 8 are, re- 
spectively, ‘Multipliction and Function Generation,” “Auxiliary Circuits 
and Computer Operation,” and ‘The Design of Complete D-c Analog- 
computer Installations.” 

The case for the use of these computers is skillfully presented and the 
authors’ concern with the accuracy of individual components is admirable. 
However, the book discusses only briefly the possibilities for large network 
computers and presumably many other types of analogs were deliberately 
excluded from the scope, for instance, a.c. analogs, continuous analogs such 
as electrolytic tanks and conducting sheets, and the many types of stress and 
strain analogs. 


F. J. M. 


1014. M. Tribus, “Intermittent heating for aircraft ice protection with 
application to propellers and jet engines,’” ASME Trans., v. 73, 1951, 
p. 1117-1130. 


The author uses an electrical analog computer for the study of deicing of 
airplane parts. Thermally speaking, the problem is that of an idealized fin in 
transient heat flow. Heat is generated in the body which is to be deiced, for 
example the propeller, by means of electric heating wires; the heat is lost to 
the surrounding environment across a boundary conductance. The paper 
deals in considerable detail with the problems of deicing but mentions only 
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very briefly the computer used. No reference is made to the large amount of 
literature on the paper. 

As far as computing technique goes the main point of interest is an elec- 
tronic circuit used to represent the boundary conductance which is not con- 
stant but rather a function of temperature. This nonlinear condition has 
heretofore on other computers been represented in discrete finite steps. The 
author shows in figure 8 a circuit for representation of such boundary con- 
ductance. The check of the computations with actual experiments is rather 
unsatisfactory possibly because of poor assumption of physical constants of 
the system; the correction of assumptions determines in analog computers 
the validity of the result. The author, in designing the computing circuit, 
disregards a number of influences (for example the thermal resistance of the 
ice). 

Victor PASCHKIS 
Heat and Mass Flow Analyzer Laboratory 
School of Engineering 
Columbia University, New York, N. Y. 


1015. F. Werner, “Further remarks on intermittent heating for aircraft 
ice protection,” ASME Trans., v. 73, 1951, p. 1131-1137. 


The paper deals with the deicing of a propeller of ovoid cross section. 
The butt end of this ovoid will be called A, the extension, B. Heat is gen- 
erated by electric heaters only in the region A and is conducted to the 
region B. A cross section through the propeller in the region A shows 
(proceeding inwards) the steel blade (extending all the way to B), a layer of 
nylon, the heater, another layer of nylon, and finally sponge rubber. Dis- 
regarding the heat flow along the propeller length, the problem is two di- 
mensional. 

The greatest problem in solving two dimensional problems on analog 
computers, one on which little information is available, is that of how to sec- 
tion the body in which heat flow occurs. The authors disregard the thermal 
resistance across the thickness of the steel, and along the length of thenylon 
and heater layers. The sponge rubber is represented by a number of parallel 
sections, ending in a fictitious center with zero volume (Fig. 3 of the paper). 
This design is not discussed or analyzed; thus, the most crucial problem, 
from a computational view-point, is not dealt with in the paper. Regarding 
deicing, the paper shows the desirability of using high rates of energy produc- 
tion, which results in a lower total heat consumption than heating at a lower 
rate. 

Victor PASscCHKIS 
Heat and Mass Flow Analyzer Laboratory 
School of Engineering 
Columbia University, New York, N. Y. 


NOTES 


143.—ANALYTICAL APPROXIMATIONS. [EDITORIAL NOTE: In Note 139, 
MTAC, v. 6, p. 251-253, Ceci. HastinGs has described the RAND Collection 
of Illustrative Approximations. The interest that these approximations have 
aroused during the past year is considerable. It is hoped to publish as Notes 
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from time to time additional examples of such approximations contributed 
by our readers. To encourage this hope, Mr. Hastings is submitting a dozen 
new examples, prepared with the assistance of Mr. JaMEs P. WonG and Mrs. 
Davin K. Haywarp. These differ from the RAND Collection in form, es- 
pecially since they do not give illustrative error curves. For convenience in 
future references we are numbering these approximations consecutively. ] 


(1) Square Root: To better than 1 part in 12 over .1 < x < 10, 
xt = (1 + 4x)/(4 + x). 

(2) Pearson Cosine Transformation : To .003 over 0 < x < 1, 
r(x) = cos (x/(1 + x#)) = (— 1 — 4x + 5x*)/(1 + 8x + 62”). 
r(x-!) = — r(x) can be used to obtain function values over 1 < x 
< 

(3) Common Logarithmic Function: To better than .005 over .1 < x 
< 1, log x = — .076 + .281 x — .238/(x + .15). 
This approximation is the result of a request for a very simple 
formula to use in the reduction of certain data. 

(4) Common Logarithmic Function: To better than .000,004 over 
10, 
log x = 4 + .86857y + .29059y* + .15783y* + .20269y’, 
where y = (x — ¥10)/(x + Vi0). 

(5) Common Logarithmic Function: To better than .000,000,015 over 
log x = 4 + .8685888y + .2895497y* + .1731159y> + .1314381y’ 
+ .0547562y® + .1832415y", 
where y = (x — V¥10)/(x + 10). 

(6) Inverse Tangent: To better than .005 over — 1 < x <1, 

arctan x = x/(1 + .28x?). 


This approximation is the result of a request for a very simple for- 
mula to use in the reduction of certain data. 


(7) Descending Exponential Function: To better than .000,000,11 
overrO<x< ow, 
= (1 + ax + + + + 
where a; = .125,000,204, a2 = .007,811,604, a; = .000,326,627, 
a, = .000,009,652 and a, = .000,000,351. 
(8) Incomplete Gamma Function Type Integral: To better than 
.000,000,1 over 0 < x < 1, 


F(x) = f 


1 
= [.219384+x(.024717+.000803x) J/[1+x(.558651 +.090584x) J]. 
F(x) + (x + 1)F(x + 1) =e can be used to obtain function 
values outside the range indicated. 
(9) Exponential Integral of Negative Argument: To better than 
.000,000,1 over 10 < x < ~, 


xe* te~‘dt 
= [1.15198 + x(4.03640 + x)]/[4.19160 + x(5.03637 + x) ]. 
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(10) Segmental Area Function : To better than .0012 over — 1 < x < 1, 
A(x) = (1 — at = 2.0083x — .4160x2 + .1604x* — .1808x". 


In terms of elementary functions, A(x) = arcsin x + x(1 — x*)!. 
(11) Segmental Area Function : To better than .00016 over — 1 < x < 1, 


A) fia ~ 


= (1.99916x — 2.39484x* + .58673x°)/(1 — 1.03472x* + .15634x*). 
(12) Segmental Area Function: To better than .000,016 over —1<x<1, 


= x(1.999872 + 4.143151y — 3.153670? — 1.430807n*)/ 
(1 + 2.901498 — 1.811287y? — 1.098016’), 
where » = x*/(5 — 4x). 
Cecit HAsTINGs, JR. 
RAND Corporation 
1700 Main Street 
Santa Monica, California 


144.—ZEROS OF THE DERIVATIVE OF BESSEL FUNCTIONS OF FRACTIONAL 
ORDER. The NBS Computation Laboratory' has published extensive tables 
of Bessel functions of fractional order, J,(x), + »v = }, 3, 4, 3, and zeros of 
J ,(x) have been tabulated by ABRAMOWITZ? and by the Computation Labora- 
tory when the latter was known as the Mathematical Tables Project.* (The 
zeros in the last two sources are also given in the first, p. 384-385.) This 
present note gives the first seven or eight positive zeros of the first derivative 
of those Bessel functions. Following standard notation, j,,, will be used to 
denote the s-th positive root of J,’(x) = 0. The zeros j;,, were obtained from 
the published values' of J,(x), and this table includes all such zeros within 
the range of tabulation of J,(x) itself (i.e., not exceeding 25). 

These zeros were computed to the maximum accuracy obtainable from 
the NBS tables. All entries are given here to seven decimal places; although 
the seventh decimal place is not absolutely guaranteed, it has a high prob- 
ability of being correct. The zeros were computed using the 5-point case of 
two different formulas for inverse interpolation for the derivative, given in 
SaLzER‘ on p. 214 and p. 215. They were also checked by (a) calculating an 
approximate expression for the error in j,,,, (b) recomputing the last 
jas given here by the asymptotic formulas for each v (see below), and (c) 
computing the seventh divided difference of 7,, as a function of », using a 
formula in Salzer.§ (This divided difference check was not fully applicable to 
the first few zeros.) 

For zeros > 25, the following asymptotic formulas, whose coefficients 
were calculated from the general expression for j,,, in WATSON,® will give at 
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least seven decimal accuracy : 


= y — .65625000y-! — .54931641y— 


__. _ f .39269908 »=— 
jhe = y — 5972222291 — .41380530y-* 
.26179939 »= M 
where y = — y=? 
4 
+v=4, = y — 43055556y— — .07507073y—* 
+v=3, je = y — .40625000y — .03108724y— 
.39269908 »=—} 
where y = + { — 196349540) 
: To obtain j,,, for any other » that is less than one in absolute value, the 
present table may be used in conjunction with a special table of interpolation 
iA coefficients,’ p. 393-413. The user is cautioned that for interpolation as well 
7 as for forming divided differences, the values of 7,,, for vy < 0 are not contin- 
ued into for > 0, but into j,, 
oi Mrs. RutH E. CaPuANo and Miss Mary M. Dun tap assisted in the 
computations. 
43 HERBERT E. SALZER 
NBSCL 
Table of 7, 
, 1 2.47861 49 2.65267 49 3.27468 22 3.41838 81 1 
i 2 5.77630 68 5.92026 00 6.47892 00 6.61491 38 2 
: 3 8.95866 64 9.09725 75 9.64204 42 9.77606 19 3 
G 4 12.11945 77 12.25581 13 12.79457 06 12.92770 62 4 
s 5 15.27226 12 15.40738 64 15.94278 34 16.07542 28 5 
6 18.42121 33 18.55556 21 19.08881 57 19.22113 82 6 
7 21.56801 08 21.70182 42 22.23359 29 22.36569 56 7 
8 24.71348 00 24.84690 21 8 
s y=i v=} v=} v=} pri 
1 1.51433 70 1.40121 80 0.90999 85 0.76906 15 1 wa 
2 4.97223 54 4.85063 49 4.35291 38 4.22515 79 2 
3 8.16610 90 8.04140 90 7.53529 41 7.40675 25 3 
+ 11.33027 35 11.20403 00 10.69360 09 10.56453 27 4 
5 14.48452 05 14.35735 04 13.84430 89 13.71489 82 5 an 
6 17.63422 27 17.50643 40 16.99164 33 16.86199 56 6 int 
7 20.78145 96 20.65322 86 20.13718 52 20.00736 48 7 
8 23.92720 84 23.79864 53 23.28166 09 23.15170 93 8 Na 
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1NBSCL, Tables of Bessel Functions of Fractional Order. V. 1, New York, Columbia 
University Press, 1 


8. 
Ps. |. eens “Zeros of certain Bessel functions of fractional order,”” MTAC, v. 1, 
p. 354 
3 Mathematical Tables Project, National Bureau of Standards, “More zeros of certain 
Bessel functions of fractional order,’’ MTAC, v. 2, p. 118-119. 


*H. E. Sauzer, “Formulas for finding the argument for which a function has a given 
derivative,” MTAC, v. 5, p. 213-215. 


5H. E. Sauzer, “The checking of functions tabulated at certain fractional points,”’ 
v. 2, p. 318-319 


N. Watson, A Treatise on the Theory of Bessel Functions. Cambridge University 
Press, 1944, p. 507. 


145.—AN EXAMPLE IN THE USE OF THE DIFFERENTIAL ANALYSER. In a 
recent article SPRAGUE! has discussed, as an example, a differential an- 
alyser setup for solving the equation 


— ww’ — wt = 0, 


4 
dt 


Figure 1 


primes denoting differentiation with respect to ¢. There are several simpler 
ways of setting up this equation. One is to take the once-integrated form 


f wd(w + + w'(0) 


and so place it on the analyser in accordance with Figure 1. This uses three 
integrators in lieu of the six shown in Figure 3 of Sprague’s paper. 

National Physical Laboratory J. G. L. MicHEL 
Teddington, Middlesex 
England 
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72 CORRIGENDA 

[EpiToriAL NOTE: Mr. SPRAGUE informs us that an additional integrator 
will be required in case the differential analyzer is of digital type, thus the 
above method would require four integrators. ] 


1R. E. SpRAGvE, “Fundamental concepts of the digital differential analyzer method of 
computation,” MTAC, v. 6, p. 41-49. 


146.—Two NEw MERSENNE PRIMES. The program described in Notes 


131(c) and 138 [MTAC, v. 6, p. 61, 204] has been continued. Two more ee 


Mersenne primes, 2° — 1 and 27%! — 1, were discovered by the SWAC on 
October 7 and 9, 1952. The time required for either of the tests is one hour. 
This makes 17 Mersenne primes, and a corresponding number of perfect 
numbers, now known. They are 2” — 1 form = 2,3,5,7, 13,17, 19, 31, 61, 89, 
107, 127, 521, 607, 1279, 2203, and 2281. 
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CORRIGENDA 
129, 1. 22 for — read + pi41) 


p. 1 

p. 132, 1. — 6 and — 18 for THompson read THOMSON 
p. 152, 1. —5 for 8 read 9 
p. 1 
p. 1 


6 
. 6, 
6 


. 6, p. 187, 1. 12 for QVAC read QUAC 
6, p. 189, 1. — 8 for ConNoLy read CONNOLLY 
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